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Abstract: 

Power corrections to exclusive processes are usually calculated using models for twist-four 
distribution amplitudes (DA) which are based on the leading-order terms in the conformal 
expansion. In this work we develop a different approach which does not rely on conformal 
symmetry but is based instead on renormalon analysis. This way we obtain an upper bound 
for the contributions of higher conformal spin operators, which translates into a bound 
on the end-point behavior of DA. The existence of such a bound is important for proving 
factorization theorems. For the two-particle twist-four DA we find in the renormalon 
model that the conformal expansion converges but it does not converge uniformly near 
the end points. This means that power corrections to observables which are particularly 
sensitive to the region where one valence quark is soft, may be underestimated when using 
the first few terms in the conformal expansion. The renormalon models of twist-four DA 
of the pion and the p meson are constructed and can be used as a viable alternative to 
existing models. 
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1 Introduction 



The relevant non-perturbative degrees of freedom in hard exclusive processes are described 
by hadron distribution amplitudes (DA) Ej that detail the momentum-fraction distri- 
butions of partons in the infinite momentum frame by integrating out the transverse 
momentum dependence. The leading-twist parton distributions appear in the QCD de- 
scription of hard exclusive reactions to the leading power accuracy and refer to parton 
configurations with the minimal number of constituents. The higher-twist distributions, 
in turn, are more numerous and are used to take into account a variety of effects due 
to parton virtuality, transverse momentum, and contributions of higher Fock states that 
are relevant for the description of power-suppressed corrections in the hard momentum. 
It should be mentioned that application of QCD factorization to exclusive processes be- 
yond the leading-twist approximation presents a serious challenge because of end-point 
divergences related to the contributions of soft partons. Up to now such applications 
have been mostly in connection with the so-called light-cone sum rules (LCSR) in 
which case the end-point divergences are removed by construction. More generally, end- 
point contributions have to be added and there is a hope that the precise separation of 
hard (twist-expandable) and soft (end-point) contributions can be achieved within the 
soft-coUinear effective theory, see e.g. Having in mind the existing and potential ap- 
plications it is worthwhile to have a fresh look at the theoretical description of higher-twist 
DA and the related uncertainties. The present work presents a step in this direction. 

The existing theoretical framework for the description of DA is based on the conformal 
symmetry of the QCD Lagrangian, see [S] for a detailed review. The symmetry can be used 
to separate the dependence of the hadron wave function on the longitudinal momentum 
fractions and the transverse coordinates (that are later traded for the renormalization 
scale) in very much the same way as the rotational symmetry of the potential in quantum 
mechanics allows to separate the angular and the radial dependence. The orthogonal 
polynomials appearing in the expansion of distribution amplitudes PI |2] are nothing but 
the irreducible representations of the coUinear conformal group and play the same role as 
spherical harmonics in quantum mechanics, with the orbital angular momentum replaced 
by conformal spin. One motivation for using the conformal expansion is that contributions 
with different conformal spin do not mix with each other under QCD evolution to the 
leading logarithmic accuracy ^ |2] . Another rationale [Hj is that QCD equations of motion 
(EOM) only relate terms with the same conformal spin so that any exact EOM-based 
relation between DA can be satisfied order by order in the conformal expansion. It follows 
that any parametrization of DA based on a truncated conformal expansion is consistent 
with the EOM and is preserved by evolution. 

Assuming that DA are dominated by first few terms having the lowest conformal spin, 
the conformal expansion provides a practical framework for constructing models for DA 
[HIHlinillDlinillllllllHIIlllISl which are both consistent with the QCD constraints and 
involve just a few non-perturbative parameters. This approach has been used extensively 
for phenomenology. Its main limitation is that just the first one or two terms in the 
expansion are included; the increasing number of parameters at higher conformal spins 
makes this program impractical otherwise, so the assumption that the conformal expansion 
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converges fast is absolutely essential. 

From the theoretical point of view, convergence of the conformal expansion is ex- 
pected because anomalous dimensions of conformal operators are rising logarithmically 
with the spin^ and, therefore, higher-spin contributions are suppressed at asymptotically 
large scales. However, this suppression is numerically weak and not sufficient to guar- 
antee convergence at the scales of practical interest. For leading twist this question can 
eventually be decided by experiment and indeed the qualitative success of quark counting 
rules and, more importantly, the CLEO data on the 71(77)77* transition form factor 
strongly suggest that the meson DA are not very far from their asymptotic shape. For 
higher twist, using experimental data to constrain DA seems totally unrealistic and up 
to now there had been no arguments whatsoever whether truncation of the conformal 
expansion is legitimate in this case or not. 

The aim of this work is to construct alternative models for twist-four meson DA that 
do not rely on the conformal expansion and present plausible bounds for the higher-spin 
contributions to these DA. To this end, we approach the problem from a different angle, 
using the concept of renormalons J7j. Owing to the breaking of scale invariance in the 
quantum theory through the running of the coupling, operators of different twist mix 
with each other under renormalization. Independence of a physical observable on the 
factorization scale implies intricate cancellations between different twists — the so-called 
cancellation of renormalon ambiguities — and the existence of these ambiguities can be 
used to estimate power-suppressed corrections in a similar way as the logarithmic scale 
dependence is used to estimate the accuracy of fixed-order perturbative calculations. 

Most applications of renormalons have so far been in inclusive or semi-inclusive cross 
sections |17| ITHj , although some work has already been done in the context of exclusive 
processes. This includes analysis of the large-order behavior of the Brodsky-Lepage evo- 
lution kernel [1^], and studies of infrared renormalons in specific processes, e.g. the 7*77ro 
form factor [W, the pion electromagnetic form factor [21], and deeply- virtual Compton 
scattering Here we do not consider a specific process but instead develop a general 
framework for estimating higher-twist contributions to exclusive processes involving pseu- 
doscalar and vector mesons by constructing a renormalon model for DA, continuing the 
work by J. Andersen We trace the cancellation of renormalons and in this way es- 
tablish an explicit connection between previous renormalon calculations and the operator 
product expansion (OPE) for exclusive processes, for cases that the latter exists. We then 
use the renormalon model to study the convergence of the conformal expansion. 

In order to explain how the concept of renormalons becomes useful for constructing 
a consistent model for higher-twist DA let us discuss a concrete example. Consider the 
following matrix element: 

(0 |T{rf>2)4^75«(^i)}| ^+(P)>^2,1/|^.| = (1.1) 

1 



i{pA) due 



'iupxi—iupx2 







L^o 

where u = 1 — m, /tt is the pion decay constant and the fields are taken at a small 



•l-This property is proven for twist-two and twist-three operators and is probably true for all twists. 
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space-like separation: A = xi — X2 with < 0, |A^| ^ I/Aq^j^. Also, /i stands for 
the ultraviolet renormalization scale and for this discussion we assume that /i^ ~ Vl^^l 
to avoid large logarithms. The expansion in the square brackets on the right-hand side 
of Eq. dHH) is nothing but the OPE: C{v,u, A'^, /j^j,) = 6{u - v) + O(a^) is the twist- 
two coefficient function, 0^(t>;/i|,) is the standard, twist-two pion DA, and A^(7i(m; 
represents the twist-four contribution; fip is a factorization scale. The function gi{u; iJ?p) 
can be interpreted as the distribution of the transverse momentum (squared) of the quark 
in the pion |Tni E] and is one example of the higher-twist DA that we will be interested 
in. 

The conformal expansion of gi{u\ yu|n) has been constructed in [HI using the expansion 
of three-particle DA involving a quark, an antiquark and a gluon, and EOM. In Sect. |3] 
we shall explain how this is done. Here we just quote the result: 

^i(u;/i|) = |^(-^=3)(/i^)[u2u2] +g^^=^\^il)^uu{l3uu + 2) (1.2) 

+ 2 (6 + 3 M + 1) ln(u) + 2 (6 + 3 M + 1) uHn{u) + ■ ■ 

where the coefficients g^'^^fJ^'p) renormalize multiplicatively, and J = 3, 4, ... is the confor- 
mal spin. The corresponding logarithmic scale dependence cancels against the scale de- 
pendence of the twist-four coefficient function which is not shown in Eq. (jl.lj) for brevity. 
To the leading logarithmic accuracy and with the usual choice = — l/l^l ^^is 
function reduces to a constant which can be included in the definition of gi{u] fij^). The 
dots in Eq. ()1.2|) stand for the contribution of conformal spins J > 5, which were usually 
neglected. In this paper we want to study their significance. 

To understand the role of renormalons one needs to carefully examine the separation 
made in Eq. (Il.lj) between twist two and twist four. For a qualitative discussion it is 
convenient to use a hard cutoff 1/|A^| ~ ^ ^ A^ for factorization: loop momenta 
k"^ > contribute to the coefficient function C{v, u, A^, /i|,) while momenta k"^ < con- 
tribute to the DA. Upon expanding the coefficient function near the light-cone | A^|yu|. <^ 1, 
one obtains: 

C{v, u, A^, fip) = (^1 + ci^s + C2a1 H j - d{v, u)/i|.A^ ■ ■ ■ , (1.3) 

where = a;s(/i^ = 1/|A^|). It is the usual perturbative series with coefficients q = 
Cj (u, m; In (/i|,//i^)) depending logarithmically on the scales, and the d{v,u) term repre- 
sents the leading power correction that arises because the low-momentum regions are cut 
off. Since the left-hand side of Eq. p.l|) does not depend on fip, any such dependence 
in C(f , M, A^, /i|,) should cancel within the square brackets on the right-hand side. In 
particular, the logarithmic dependence of C(f , m, A^, /i|,) on is canceled by that of 
the twist-two distribution amplitude (p-^iv] 11^). The cancellation of power dependence 
must involve the twist-four term A'^ gi{u; fij^) , so it is expected that in this factorization 
prescription 

gi{u;fil) — > fil / dvd{v,u)(j)^{v;fip)+gi{u;fip), (1.4) 
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where the second term depends on at most logarithmically. Indeed, upon renormal- 
ization the relevant twist-four operators show not only logarithmic ultraviolet divergence 
which is related to their anomalous dimension, but also quadratic ultraviolet divergence. 
Using a hard cutoff, the dependence of twist-four on /ij^ is indeed that of Eq. ()1.4j) : 
the twist -four operators mix with the leading-twist such that the dependence on in 
Eq. (jLip cancels out. Clearly the function d{v,u) can be computed either by considering 
the dependence of the twist-two coefficient function on fi'jp, i.e. its sensitivity to small 
loop momenta, or by considering the dependence of twist-four operators on i.e. their 
sensitivity to large loop momenta. Cancellation of /il^ to power accuracy requires that 
these two regularizations will be done using the same prescription. 

In practice, a hard cutoff is difficult to implement. Usually, dimensional regularization 
is used instead. In this case power terms in the coefficient function (Eq. p.3p ) do not ap- 
pear. However, the coefficients Cj computed in a MS-like scheme diverge factorially with 
the order i. The factorial divergence implies that the sum of the perturbative series is only 
defined to a power accuracy and this ambiguity (renormalon ambiguity) must be compen- 
sated by adding a non-perturbative higher-twist correction. A detailed analysis shows 
fri\ that the asymptotic large-order behavior of the coefficients (the renormalons) is in 
one-to-one correspondence with the sensitivity to extreme (small or large) loop momenta 
and that infrared renormalons in the leading-twist coefficient function are compensated 
by ultraviolet renormalons in the matrix elements of twist-four operators. At the end the 
picture described above re-appears: only the details depend on the factorization method. 

Returning to Eq. ()1.4|) we observe that the quadratic term in /xp is spurious since its 
sole purpose is to cancel the similar contribution to the coefficient function. It does not 
contribute, therefore, to any physical observable. The idea of the renormalon model [TTllTSj 
is that, with a replacement of fip hj a suitable non-perturbative scale, this contribution 
should be of the same order and have roughly the same functional form as the physical 
second contribution on the right-hand side of Eq. ()1.4|) . which is the only one of interest. 
Assuming this "ultraviolet dominance" fI^ |23 EI we get the following model: 



dv d{v,u)(f)T,{v; fxj^), (1.5) 



9iiu;fil) ~ fcA^cD / 

^0 

where by explicit calculation (see Sect. one finds for d{v, u) at leading order 

div, u) = — M + (f — m) In f 1 \ 9(v > u) 

L V t>/ J 

1 — u + iu — v)\n [ 1 



1 

+ 



{l-vf 



1 -u 



l-v 



e{v<u). (1.6) 



The overall coefficient in this model k = 0{1) can be fixed by the normalization inte- 
gral Jq du gi{u; fi"^) corresponding to the matrix element of a local operator which can 
be estimated by QCD sum rules [TUl E] or calculated on the lattice. To the accuracy of 
Eq. p.6|) the logarithmic //j? dependence on the left-hand side and on the right-hand side 
of Eq. (jl.5p do not match: one expects this model to be relevant at low scales of the order 
of a few times Aqcd- 
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The "ultraviolet dominance" assumption used to derive Eq. (jl.Sp , is in fact sufficient to 
derive the full set of two- and three-particle DA of twist-four in terms of the leading-twist 
DA. It is important that this approximation is fully consistent with the OPE and respects 
all constraints imposed by EOM. On the other hand, it does not assume any hierarchy 
of contributions of the increasing conformal spin and indeed the expressions in Eq. (jl.2|) 
and Eq. ()1.5|) are very different for any reasonable choice of the leading-twist DA. Since 
twist-four anomalous dimensions are not taken into account, one should expect that this 
model overestimates contributions of conformal operators with high spin. It is therefore 
complementary to the usual models based on using a ffist few terms with the lowest spin 
and provides an upper-bound estimate of the neglected contributions. 

The presentation is organized as follows. In Sect. 2 we formulate our program in more 
precise terms, introducing the relevant techniques (Borel transform) and the systematic 
approximation (large-iVj expansion) that will be used throughout the rest of the work. 
The expression in Eq. is derived and the cancellation of the renormalon ambiguity 

for the matrix element of the type in Eq. (jl.lj) is demonstrated by an explicit calculation. 
The renormalon model of the pion DA of twist-four is presented in Sect. 3 and in Sect. 4 
we discuss its conformal expansion. The generalization of these results for the case of 
vector mesons is considered in Sect. 5 while Sect. 6 is reserved for the conclusions. 



2 Cancellation of renormalon ambiguities in the OPE 
2.1 Definitions 

In this section we compute the renormalon ambiguity of the leading-twist coefficient func- 
tion in the simplest exclusive amplitude involving a pseudoscalar meson, and demonstrate 
how the unique result is restored by adding the twist-four contributions in the OPE. To 
begin with, let us set up the necessary definitions. 

We choose to consider the gauge-invariant time-ordered product of two quark "cur- 
rents" at a small (but non- vanishing) light-cone separation, which can be parametrized in 
terms of two Lorentz-invariant amplitudes Gi{u, A^; /i^) and G2{u, A^; /i^) defined as 

{0\T{d{x2hu'y5[x2,Xi]u{xi)}\n+{p))^2 = (2.1) 







Gi{u, A'; /i')p, + G2{u, A'; fi') { ^^A, - 



Here A = X1—X2 with |A^| ^ A^ < playing the role of the hard scale, u = 1 — u 

and /i^ is the ultraviolet renormalization scale. We use the notation [x2, xi] for the Wilson 
line connecting the points X2 and xi, 



[x2, xi] = Pexp 



-ig / dtA^A^'{x2 + tA) 



(2.2) 



Note that the /i^ dependence comes entirely from the wave-function renormalization of the 
quark and the antiquark fields in Eq. ()2.1|) and can be removed by adding the corresponding 
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Z-factors. Up to this additional renormalization, the functions G'i_2(m, A^; /x^) can be 
viewed as physical amplitudes. Their advantage is that they are simpler than exclusive 
amplitudes which are relevant for phenomenology (such as the two-photon one) and at 
the same time retain all the features that are important in the present context. 

The asymptotic behavior of Gi{u, A^; n"^), i = 1,2 in the light-cone limit A^ — > 
with A • p fixed can be studied by the OPE, schematically 

G.iu, A') = Cf ® + E ^^'^ ® + ^(^')' (2.3) 

where Cf^ are the coefficient functions and 0^*-' are the pion DA given by vacuum-to-pion 
matrix elements of renormalized non-local operators, t refers to twist and the summation 
goes over all independent contributions of given twist. The convolution is defined as 

®0 = / dvCi{v,u,A^,fi^,fil)(f){v;fil), (2.4) 

where u, v have the meaning of the light-cone momentum fractions and we have indicated 
the dependence on the factorization scale fijp and other variables. The Lorentz structures 
in Eq. ()2.H) are chosen such that the coefficient functions Cj depend on A^ logarithmically 
as follows from power counting. 

The leading-twist pion DA 0,r(w) = (j)^'^\u) is defined as usual by 

{Q\d{-z) ^^,[-z,z]u{z)\iT+{p)) =iUpz) [' due-''^'^^-^Un{u-fil), (2.5) 

^ Jo 

where the normalization condition is 

du (j)T,{u] fj^p) = 1. (2.6) 







Here and below 2;^ is the light-like vector, z"^ = 0. Although we will usually be using 
covariant notation, it is sometimes convenient to refer to light-cone coordinates where 
p = (p+, 0, 0) and x = (0, x_, x±). The light-cone limit corresponds to x± — > with x_ 
fixed such that z = (0,x_,0). 

In the present paper we will consider the leading-twist coefficient function to all or- 
ders in a<j(/x^) but restrict ourselves to the leading order in higher-twist contributions: 
CW(u, v) = 6{u -v) + 0{a,). To this accuracy there is a single twist-four DA contribut- 
ing to each of the functions 6*1.2: 



G,{u, A^) 
G2(u, A^l 



1 + cias + C2a1 + 



cias + C2C?^ + 



b^\A^f{u-ix^p)\--- 
A^f{u-^l)^---. (2.7) 



In physical terms ^[^''(m; [i^p) and <^^^{u] /i|^) correspond to contributions of valence-quark 
transverse momenta and of the "wrong" components of the quark spinors, respectively. In 
the notations of Ref. 0| 

9\{u) = (l)^i\u) and 5(2 (m) = ^0^^^(m). 
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For an operator definition, it is convenient to use the concept of non-local operator OjO^the 
light cone [2^1 which is defined as the generating function (formal Taylor expansion) for 
renormalized local operators. For example, 



d{—x) ^■y5[—x,x]u{x) 



E 



2fc+i 



d'ju D 



Dui u — Traces 



twist— 2 



I rrp 1 twist— 4 I rrp n twist— 6 , 

-|-[iracesj^2^ + [iracesj,,2 + 



(2i 



where Dui=Diy^ — Dux is the covariant derivative, see [211 I2Z| for details. Note that 
the twist expansion of a non-local operator corresponds to its expansion over irreducible 
representations of the Lorentz group. This involves rearrangement of traces and sym- 
metrization over groups of indices if necessary, c/. j2H] • This compact notation is widely 
used in UHl HB 123 CSl CSl • In particular, we have 



(0| d{x2)'^ulb[x2,xi]u{xi)} 



TT 



\P)) 



-I Mf 



'o 



?)(«;/i^) + 0(A2 



(2.9) 



Despite the similar appearance, beyond the tree level Eq. (j2.1|) and Eq. (j2.9p describe 
different objects. As follows from the Taylor expansion Eq. (j2.8p the non-local operator 
on the left-hand side of Eq. ()2.9|) is an analytic function of A^ at A^ — > 0. It corresponds 
to the analytic part of the amplitude in Eq. ()2.1|) P7] while the coefficient functions (by 
definition) take into account singular contributions. The most striking difference is that 
the non-local operator in Eq. (j2.9|) does not include the whole twist-two part of (^2(^5 ^^) 
in Eq. (j2.ip . which is ~ 1/A^. This contribution was overlooked in [221 • The non-local 
operator also does not include any of the coefficient functions appearing in Eq. ()2.1|) . which 
absorb all logarithms ~ ln(— A^/i|,). 

Using the EOM the two-particle pion DA (t)t\u\ /i|,) and 02^''(m; /il^) can be expressed 
in terms of the Fock components involving an extra gluon field. Following ^ we define 
the three-particle pion DA of twist four 

(0 \d{-z)[-z,vz\-iy'^^gG^p{vz)[vz,z\u{z)\Ti''{p)) = 

= Ujva, e-^^'^(-i-"2+-^^)|^(p^;2,-Pp;2^)<l>ll(ai,«2,«3) 



+ 



Pp dp^v 



Zf^Pv 

pz 



Ppi 9pv 



ZpPv 

pz 



<l>_L(ai,a2,a3) 



(2.10) 



where the longitudinal momentum fraction of the gluon is and the integration measure 
is defined as 



j T>ai = J daida2da^ 5(1 — ai — 0:2 — as) 



(2.11) 
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One obtains [H] (see also Appendix A in P^): 



#(n) + #(n) 



dv / dai 



l-v 



dai 



l-u 



da2 — r2$_L - («!, 02, as) 
as 



aa2 2 - $||J (ai, a2, as) 



(2.12) 



where as = 1 — ai — a2. The derivation of these relations relies on exact operator identities 
[2n] which relate integrals over v of the quark-gluon-antiquark operator in Eq. ()2.10|) to 
derivatives of the quark-antiquark operator appearing in Eq. ()2.9p . 

For completeness, we present here the definitions of the other three-particle twist- 
four pion DA. These distributions do not contribute to Eq. (j2.H) to leading order because 
of our specific choice of the (simple) correlation function, but are important for other 
applications. We will use these definitions in Sect. 3 where we construct the renormalon 
model. 

First of all, there exist another pair of DA that correspond to the substitution^ of 'j^G^p 
by iG^'P = |e'^^°^G„/3 in Eq. fTT^ . To the twist-four accuracy 

{0\d{-z) [-Z, vz]'jJgG^p{vz) [vz, z]u{z)\7r^ (p)) = 



/ Dai e 



-^Ma,~a,+asv)^^ (p^;2, - p,;^^) V]/ 1| («! , «2, as) 



+ 



Pp 9tiu 



pz 



ZpPu 

pz 



\l/^(ai,a2,as) 



(2.13) 



The functions ±,\\ and are not independent. In particular and $_l can be obtained 
as the symmetric and the antisymmetric part, respectively, of a more general DA 



(0|(i(-z) [-Z, vz]-f_j^j+igGpp{vz) [vz, z]u{z)\tt^ (p)) 



(2.14) 



Pp 9pv 



ZpPv 

pz 



Pp 9pv 



ZpPv 

pz 



, a2,asj . 



One obtains [12] 



^±(ai, a2, as) = -- [H^\ai, a2, as) + -^^^"2, ai, as)] 



$±(ai, aa, as) = -- [H^^{ai, aa, as) - i^-^^(a2, ai, as 

see also Appendix A in [5U] . 

In addition, we introduce a new three-particle DA ^^^(aj) as 

(O \d{-z)-^p-^r^gD"'Gap{vz)[vz,z\u{z)\'K''{p)) = 



ifnPpPfB / l^ai e 



-ipz{ai—a2+va3) 



^(ai, a2, as) 



(2.15) 



(2.16) 



§We are using the sign conventions for the e^''"^ tensor and the 75 matrix from Bjorken and Drell j29| . 
In particular, Tr(7''7''7"7''75) = Aie'^''"^ and CTq,^ = idap^z- 
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The Lorentz structure Pfj,Pf3 is the only one relevant at twist four. Thanks to the equation of 
motion D"^G^p = —g Qt^l/sfl where the summation goes over all light flavors, ^^^(aj) can 
be viewed as describing either a quark-antiquark-gluon or a specific four-quark component 
of the pion: with the quark-ant iquark pair in a color-octet state and at the same space- 
time point. This DA was not considered previously because its conformal expansion starts 
with a higher spin J = 5 (see below) whereas in [HI [H] only the terms with J = 3 and 
J = 4 were included. 

There exist further twist-four four-quark operators where all quark light-cone coordi- 
nates are separated, and also twist-four operators which include two gluons in addition to 
the quark-ant iquark pair. They give rise to four-particle DA and will not be considered in 
this work. Although such operators contribute to exclusive amplitudes at twist four, they 
do not appear to leading order in the flavor expansion (c/. [SI]) and can be systematically 
neglected to our accuracy. 

2.2 Infrared renormalons in the twist— two coefficient functions 

The scale separation made in Eq. ()2.3|) is arbitrary in two respects. To a given logarithmic 
accuracy the separation between coefficient functions and operator matrix elements is 
ambiguous. However, the logarithmic dependence on the factorization scale cancels 
between the coefficient function (7*^*^^^ and the corresponding DA 0^*"^) leaving their 
product invariant. This can be phrased through perturbative evolution equations. To 
power accuracy the separation between terms of different twist is ambiguous and the 
scale-independence of "structure functions" G is only restored in the sum of all twists. 

As discussed in the introduction the intuitive way to see this is to imagine using a 
cutoff scale fi'jp to implement factorization: this scale serves as an infrared cutoff for 
the coefficient functions while it acts as an ultraviolet cutoff for the DA. The twist-two 
contribution would then depend on the scale as ~ jj?p/\^(j)'^'^\ To keep G invariant this 
dependence should cancel against a term of the form j^'jp(f)'''^'' in ^j-^"*. In other words, the 
twist -four operator has a quadratic ultraviolet divergence through which it mixes with 
the twist-two operator. In dimensional regularization the power-like cutoff dependence 
does not occur, but the ambiguity still persists because the perturbative series develops a 
factorial large-order behaviour. The perturbation theory, therefore, diverges (renormalon 
divergence) and its sum is only defined to power-like accuracy. This is the renormalon 
ambiguity which we are going to address now. 

2.2.1 An all— order calculation 

In order to evaluate the infrared-renormalon ambiguity in the twist-two coefficient func- 
tions 

C[^\v, u, A^, fi"^, fi%) = S{u - v) + ciasifi^) + C2a1{fi'^) H 

(i;,n,A2,/i2,/i|) =cia,ifi^) + C2alifi^) + --- (2.17) 

we need to calculate them to all orders in the coupling. Of course, a full all-order calcu- 
lation cannot be performed. Instead, as in other applications ^7], we restrict ourselves to 
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the perturbative series generated by the running-couphng effects in one-loop diagrams, 
i.e. using QCD couphng at the scale of the gluon virtuality^. A convenient tool to perform 
this calculation is the Borel representation of the running coupling, 



Jo v-^^ 

where /3o is the leading-order coefficient of the /3 function, 

P.-'^,C.-In,. (2.19) 

A = Aq^qj-) and the exponential factor et"' originates from the renormalization of the 
fermion loop in the MS scheme. This way the calculation reduces to one loop with a 
modified gluon propagator 

— > -^—-^ . (2.20) 

and the result for the coefficient functions takes the form of a Borel integral 

POO 

Jo 

POO 

C2{v,u,^\^i\^il)= / f/w7ei"'52(«;;T;,M,AV,/iF//"') (-A'A^)". (2.21) 

Jo 

The perturbative expansion of Ci can be recovered order by order from the expansion of 
their Borel transforms near = observing that for one-loop coupling 



(In/iVA 



2\n+l 



The fixed-sign factorial behavior of the coefficients in this expansion in high orders man- 
ifests itself in that the functions Bi develop singularities at finite values of w on the 
positive real axis (renormalon singularities) rending the integrals in Eq. ()2.21|) ill-defined. 
The imaginary part that arises when the integration contour is moved above (or below) of 
the nearest singularity (to the origin) can be taken as a measure of the ambiguity of the 
summation of the perturbative series, see |^ for details. 

The relevant diagrams are shown in Fig. ^ The calculation is presented in detail in 
Appendix^ We consider quarks on-shell and set the space-time dimension to d = 4 from 
the beginning since w 7^ acts as a regulator. The result for the (not yet renormalized) 
coefficient functions reads 

BM.u^v) = ^^^,^-{e{u > .)i [(1 +«;)/r (1 - (1 - ^ 



+ e{u <v)- \{i + w)&^ (i- -)- wf^'> f 1 - - 

V I \ V / \ V 



^This approximation is sometimes referred to as the large-/3o limit and can formally be defined con- 
sidering first the large Nf limit with UsNf is fixed, and replacing Nf by — 6/3o to recover the non-Abelian 
contributions. 



11 



X2 



X2 



Figure 1: One-loop diagrams appearing in the non-local matrix element of Eq. H2.1|) . The upper 
double line represents a path-ordered exponential (a Wilson line) connecting xi and X2, while 
the dashed line represents a dressed gluon (dressing is equivalent to summing over any number 
of fermion-loops insertions in the large- A^y limit). 



1 — w 
1 + w 

1 

l-2w 



and 



B2{w;u,v) 



CfT(I-w) 



f3o T{2 + w) 
where we used a short -hand notation 



(n > v) 



5{u — v) 



u 

[u> V] \- 



u 



w+1 



+ d{u < v) 



+ ti(v > u) [- 



l + w 



In particular 



/(-=o)(/5) 



1-/3 



/(-=i)(/5) = l-/3 + /51n/?. 
The "+" prescription is defined as usual by: 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



The terms ~ /^"'•' in Eq. (j2.22j) correspond, in Feynman gauge, to the contribution of the 
vertex correction (Diagram 1) in Fig. ^ and its symmetric counterpart, the contribution 



12 



in the third hne in Eq. ()2.22|) and the entire Eq. ()2.23p originate from the box diagram 
(Diagram 2) and the remaining term ~ 1/(1 — 2w) in the last hne of Eq. ()2.22j) stands for 
the self-energy insertion in the Wilson line (Diagram 3). 



2.2.2 Singularities of the Borel transform 

The answer for Bi{w,u,v) in Eq. ()2.22|1 has a simple pole at w = 0. This singularity is 
expected [T7j and has to be removed by the subtraction of ultraviolet (UV) logarithmic 
divergences due to the wave-function renormalization of the quark fields, and of infrared 
(IR) logarithmic divergences that correspond to the renormalization of the leading-twist 
DA. To this end, note that vanishing of the Diagram 4 in Fig. ^ is a result of an exact 
cancellation between IR and UV divergent integrals. Upon introducing a scale which 
regulates one singularity, the other one will appear as a 1/w pole. Schematically 





'1 " 




"1" 


( 


w 


UV 


W 



[2.27) 



In addition to renormalizing Diagram 4, ultraviolet subtraction removes the 1/w pole of 
Diagram 3 in Fig. ^ 

Having performed the UV renormalization the remaining 1/w singularities from Dia- 
grams 1, 2 and 4 are removed upon performing IR factorization. The counter-term which 
is to be subtracted from the DA (p^"^^ and added to the coefficient function Bi{w;u,v) is 

^.(.,; ^ II + (2.28) 

2/?o u; 

where K{u, v) is obtained as the limit at w — >■ of the expression in the curly brackets in 
Eq. and adding the IR-pole contribution in Eq. (ITTfll : 



K{u, v) 



+ e{u >v) + ^-^ (i + ^— ) e{u < v) 

u \ u — V I 1 — n\ V — u 



(2.29) 



+ 



K{u,v) can easily be identified as the leading-order Brodsky-Lepage kernel ^112] control- 
ling the /i|.-evolution of the pion DA 0,r(M; A^f): 



(iln/i|, 27r 



du K{u,v) (j)T,{u; j2p) . (2.30) 







The terms 0{w) in Eq. ()2.28|) correspond to scheme-dependent higher-order contributions 
to the Brodsky-Lepage evolution kernel. In MS-like schemes the kernel has an expansion 
in as with a finite radius of convergence, see ^H] for explicit expressions in the large-/?o 
limit. In other words in such schemes the infrared counter-term is free of Borel singularities. 

To summarize, the subtraction of logarithmic UV and IR divergences removes the 1/w 
singularity in the Borel transformed coefficient functions in much the same way as poles 
1 / {d—4) are subtracted in renormalized amplitudes when using dimensional regularization. 
For the subsequent discussion it is important that in MS-like schemes the subtracted 
terms are analytic functions of the Borel variable w and do not infiuence the structure of 
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singularities of -Bj at w > which we are going to address now. For this reason we can 
work with non-subtracted amphtudes in what follows. 

First, we note the presence of a Borel singularity at w = 1/2 in the last term in 
Eq. (j2.22j) which comes from the self-energy insertion in the Wilson line and reflects a 
linear divergence in the UV region (UV renormalon). Such singularities are well-known in 
the context of the heavy-quark effective theory |32] in which case they reflect ambiguities 
in the non-perturbative definition of the heavy-quark mass [HHl In our case, this 
singularity is an artifact of choosing an oversimplified "exclusive process" in Eq. ()2.1|) 
where a dynamical quark propagating between the points xi and X2 is replaced by a path- 
ordered exponential. It has nothing to do with the twist expansion and will not appear 
in realistic physics applications. Therefore we will not consider this singularity further in 
this paper. 

The remaining singularities at positive integer wq = 1,2,... have IR origin and are 
called IR renormalons. They obstruct the Borel integration in Eq. (j2.21|) and render 
the sum of perturbation theory ambiguous to power accuracy ~ (A^A^)*"". Hence we 
concentrate on the IR renormalon with wq = 1 which is the closest one to the origin and 
the only one relevant for the calculation to the twist-four accuracy ~ A^. For definiteness, 
we choose the residue at the w = 1 singularity times vr = 3.14. . . as a measure of the 
ambiguity of the Borel integral (j2.2H) . i.e. 



5t, 



w 



-nf{w = 1) 



(2.31) 



Using Eq. ()2.22p and Eq. ()2.23|) we obtain the ambiguity of the twist -two approximation 
for the amplitudes Gi and G2 



5:k {Ci ® J = -cA^ A^ / dv Mv) - (i--)e{v>u) +^ /(I) (i--)e{v< 



u 



u + (v — u)\n [1 



u 



6{v > u) 



V 



u + [u — v)hi\l — — 



u 



[v < u] 



^iR {C2 ® J = cA^ A^ dv (l)^{v) I ' d{v >u)+ ' d{v < u) 
respectively, where we used Eq. fl2.25p and where the overall normalization 



(2.32) 



(2.33) 



corresponds to the convention in Eq. ()2.3H) . The given number is for Np = 3. We are going 
to demonstrate that this ambiguity is exactly canceled by the UV renormalon ambiguity 
in the twist-four DA, which is reminiscent of quadratic UV divergence of the contributing 
operators. 
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vz z 



z vz 



z 




Figure 2: One-loop diagrams appearing in the non-local matrix element of quark-antiquark- 
gluon light-cone operators. The upper double line represents a path-ordered exponential (a 
Wilson line) along the light cone connecting between —z and vz and between vz and z while the 
dashed line represents a dressed gluon (dressing is equivalent to summing over any number of 
fermion-loops insertions in the large- A'^^ limit). 



2.3 UV renormalons in higher— twist operators and cancellation 
of ambiguities 

In order to reveal the UV-renormalon divergence in the twist-four DA we consider the 
perturbative series generated by running-coupling effects in the matrix element of a generic 
quark-ant iquark-gluon operator 



d{—z)\—z^ vz)^gGaji{yz)\vz^ z\u{z) 



(2.34) 



sandwiched between quarks states with momenta q\ and g2 with ^ 0. Here V stands 
for an arbitrary Dirac structure and z"^ = 0. The relevant diagrams are shown in Fig. |21 
One obtains " 



[q2 \d{-z)TgG^fs{vz)u{z)\q,) = e-^(^^+«^)^ {g,^g^^ - g.^g^) 



l3o 



2\w 



X dg,^^^pTug, /^^(g2, (1 + v)z) - dg,T^p^^Ug, I (q^ , (1 - v)z)^ , (2.35) 

where dq^ and Uq-^ are quark spinors and I^P{q, z) stands for the momentum integral 

k^{q + k)P 



d^k 



-ikz 



(27r)4(fc2)(i+-) ^q^ky 



(2.36) 



which is UV divergent at w = 1. Performing the integral and extracting the w = 1 pole, 
we obtain: 



—I 



'"'=1 327r2ri-w) 



1 



iq^zP{l - a) - iqPz^a + g^P + -a(l - a)q^z^zP 



(2.37) 



"in the following for brevity we do not show the Wilson lines in the non-local operators. 
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As a consequence of the w = 1 singularity the Borel integral in Eq. ()2.35p is ill-defined. 
Using Eq. ()2.3ip to quantify the ambiguity and specifying for the relevant Dirac and 
Lorentz structures the final result can be brought to a form of operator relations (c/. 

CUES!) 



{z^ d{-z)Yl59Gf,u{vz)u{z)} =2icA^ / da{l-a) d{-y)^'^^u{z) - d{-z)^-i^,u 



{z^ d{-z)^j5gG^p{vz)u{z)} = -2icl\^ Zp 



da a 



z^d{—z)Yi9Gpi,{yz)u{z 



5uv i z^d{-z)^igGp,p{vz)u{z) 



2zcA^ / dail + a) 



d{,-y)h^u{z) - d{-z)^-fr^u{y) 
d{.-y)hbu{z) + d{-z)^-i^u{y 



2i chrzr, 



da a 



d{-y)h^u{z) + rf(-2;)/75M( 



(2.38) 



where y = z{l — (1 + v){l — a)) and y = z{l — (1 — v){l — a)), and c is the constant 
defined in Eq. ()2.33|) . In order to arrive at these expressions we performed integration by 
parts over a in order to remove factors of g ■ z and then converted the results into operator 
relations in configuration space. The relations in Eq. (j2.38j) can be viewed as the mixing 
under renormalization between the twist-four quark-gluon-antiquark operators and the 
twist -two quark-antiquark operator and the first two of them were derived in [22] • In a 
similar manner — see Appendix^ — we obtain one more operator relation 



5w {d{-z)z''gD''G^f,{vz)h5uiz)} 



-AicA^ 



da a- 



d 

da 



d{-y)i-f^u{z) + -^d{-z)i-i^u{ 



l-v 



(2.39) 



Taking matrix elements of Eqs. ()2.38j) and ()2.39|) between the vacuum and the pion state 
we derive UV ambiguities of three-particle pion DA defined in Sect. 2.1 in terms of the 
leading-twist pion DA: 



5uv{'^'±(ai,a2,a3)} = cA^ 
5uv {$||(ai,a2,a3)} = 2cA 
^uv{^±(ai,a2,a3)} = cA^ 
(5uv 1^11(01,02,03)} = -2cA^ 
(5uv {2^(01,02, 03)} = -4cA^ 



07r(Ol) _ 07r(O2) 
1 — «! 1 — 02 

O207r(Ol) _ Oi0^(a2) 

.(l-oi)2 {l-a^f 



07r(Ol) ^ 0,r(O2) 
1 — «! 1 — ^2 



O207r(Ol) Oi0^(a2) 
O2 07r(Ol) Oi 0^(02) 



1 — «! 



1 - O2 



(2.40) 
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where we used the symmetry property 07r(a) = (f)Tr{l — a) to arrive at the given expressions. 

Last but not least, we use EOM in Eq. ()2.12|) to calculate the UV-renormalon ambiguity 
in the two-particle twist-four pion DA (p^^l and find that it coincides identically with the 
IR renormalon ambiguity of the twist-two result in Eq. ()2.H2j) but has opposite sign. It 
follows that the "structure functions" Gi are unambiguous to the twist-four accuracy 

S{Gi} = 5:r, {Ci ® J + A^^uv {(p? } = , 

5 {G2} = 5:r, {G2 ® J + AXv {0^^ } = , (2.41) 

as expected. For the cancellation to hold, it is important that both the leading-twist 
coefficient functions and the matrix elements of higher-twist operators are calculated using 
the same regularization and renormalization prescription. This can be seen as a consistency 
check for the OPE analysis. 



3 Renormalon model for twist— four DA of the pion 

The UV-renormalon ambiguities in the twist-four DA should be viewed as indicative of 
the size and the momentum-fraction dependence of "genuine" non-perturbative effects. 
We define the renormalon model for twist-four DA of the pion by taking the functional 
form of the corresponding UV-renormalon ambiguities, replacing the overall normalization 
constant cA^ by a suitable non-perturbative parameter. The crucial observation is that al- 
though the absolute normalization of the renormalon ambiguity in Eq. ()2.31|) is essentially 
ad hoc, the relative normalization for the different DA in Eq. ()2.40|) is meaningful since the 
running-coupling calculation satisfies all constraints imposed by Lorentz symmetry and 
the EOM. Therefore, the renormalon model for the entire set of twist-four DA of the pion 
has just one free parameter: the overall normalization. This parameter can be related to 
the matrix element of the local operator 



{0\d'yJgG^pu\n^{p)) 



6^ ~ 0.2 GeV^ 



(3.1) 



where the number comes from QCD sum rules The UV-renormalon ambiguity of the 
same matrix element is, on the other hand 



5uv i {0\d-fJgG^pu\7r+{p)) \ = 2/^cA [ppg^^ -p^gpv] 



(3.2) 



so in Eq. ()2.4U|) we replace 

6 

and end up with the set of three-particle DA 



cA" ^ ~ (180 MeV)' 



$_L(ai,tt2,a3) 



$ll(a;i, 0:2, 0:3) 



52 



07r(ai) _ 07r(a2) 
1 — «! 1 — 0^2 



\2 V 



a207r(ttl) ai07r(a2) 



;i -ai) 



:i-a2) 



(3.3) 
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^±("1,02, "a) 



^11(01,02,03) 



6 



07r(ai) , 07r(a2 



1 + 1 

1 — Ol 1 — 02 

52 rO20^r(Oi)_ O10^(O2) 
'3 L(l-Oi)2 + (1-02)2 

25^ ro2 0^(oi) oi 0^(02) 



1 — Oi 



1 - 02 



(3.4) 



Note that we have made the same replacement in all five DA. It must be so because these 
DA are all related — see Appendix El Using the EOM, Eq. (Q, (c/. Eq. ^^), 



2 /•! 



8_ 
6 Jo 



u + (f — u) In ( 1 — 



u 



+ 



M + (u — 'L>) In ( 1 



6{v > u) 

e{v<u)], 



6{v < u) 



(3.5) 



which completes the calculation. 

The given expressions are valid for an arbitrary leading-twist pion DA For 
practical applications it may be worthwhile to choose the asymptotic expression 



4>n{u) = 6u{l — u) 

which is known to provide one with a reasonable accuracy ^Hl- With this choice 
<1>_l(oi, 02, 03) = 5'^[ai - 02] , 

$11(01,02,03) = 25^0102 



(3.6) 



1 — Oi 1 — O2 



\E'_l(oi, 02, 03) = (5^ [oi + 02 
\E'lj(oi, 02, 03) = — 25^0102 



1 1 

+ 



1 — Oi 1 — O2 



V(Oi, O2, O3J 



0, 



(3.7) 



Note that 2,^(01,02,03) is vanishing in this approximation, but in general it is not zero 
and has to be taken into account if corrections to the asymptotic pion DA are included. 
For the two particle DA one gets 



f)(M) = 5''{u ln(M) -Li2(M) 



+ u 



\n{u) - lA2{u) 



uu + 



-I 
6 /' 



{u) = 52 ln(M) + \Yi{u) + mm} , 



(3.^ 



where Lia(x) = ^J^ia:"/n". In the next section we shall compare these expressions with 
the model of [HI El based on conformal expansion. 
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4 Conformal expansion 



4.1 General formalism 

For fields "living" on the light-cone the conformal transformations reduce to the three- 
parameter group SL[2, M.) with the algebra of hyperbolic rotations, see ^ for a review. 
The conformal transformations for quantum fields are governed by their conformal spin 
which is defined as 

J = lis + i), (4.1) 

where i is the scaling dimension, £ = 3/2 for quarks and £ = 1 for gluons, and s is the spin 
projection on the light cone. For Dirac spinors (quarks) the different spin components can 
be separated with the help of projection operators q± = Il±q where 

n+ = ^7-7+, n^ = ^7+7-, n+ + n_ = i. (4.2) 

The g+ and g_ fields have spin projections s = 1/2 and s = —1/2, and hence different 
conformal spins j = 1 and j = 1/2, respectively. Similarly, the gluon field strength has to 
be decomposed in different spin components: G+_l has spin projection s = 1 and conformal 

spin j = 3/2; G±± and both have s = and j = 1; finally G-± has s = —1 and 

J = 1/2. 

Conformal expansion of DA presents an example of the classical problem of spin sum- 
mation in quantum mechanics, with the only difference being that the total conformal spin 
J of the multi-parton system is always larger or equal to the sum of spins of constituents: 

J = Jmin + N , Jmin = jl + ■ ■ ■ + jk ■ (4.3) 

The integer number N can be identified with the total number of covariant derivatives 
in the corresponding local operator. Adding a derivative increases the conformal spin by 
one unit and does not change the twist of the operator, defined as dimension minus spin. 
The generic conformal expansion for two-particle DA has the form 



Hu) = -01(2^'! - cjP^--'^--^{2u - 1) , (4.4) 

where P^'^^ [x\ are Jacobi polynomials jj35j , u and 1 — u stand for the momentum fractions 
of the parton with spin ji and j2, respectively, and the coefficients cj correspond to the 
contribution of the total conformal spin J = ji + j2 + N. The factor in front of the sum 
is called the asymptotic distribution amplitude. 

For three partons, a generic DA can be written as a double sum 



^ r(2ji + 2j2 + 2j, 

where ^ EHl 



N 



^ r(2j0r(2j2)r(2j3) ' ' ' hi'o 



Fjf ^(a.) = (1 - as)---- if5i^^-^^[l - 2^3] i^Sf 



a2 — «! 

1 - "3 



(4.6) 
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are the basis functions j3] corresponding to the total conformal spin J = ji+j2+j3+N and 
the fixed conformal spin of the (1, 2)-parton pair j = ji + j2 + The functions Yjj^'*^(ai) 
form a complete basis and are mutually orthogonal with respect to the conformal scalar 
product: 

Va, al'^^'al^'-'af-' Yff^\a,)Y^}f\a,) = Uj, 5jr 8,,. , (4.7) 



where 
AO. 



r(j+ji-j2)r(j-ji+j2) r(j-j+j3)r(J+j-j3) 



r(j-ji-j2+i)r(j+ji+i2-i)(2j-i) r(j-j-j3+i)r(j+j+j3-i)(2J-i)" 

(4.^ 



Using this orthogonality property the coefficients in the expansion of any DA can be 
obtained by projection: 



^ m,m2n)T(2n) 1 /■■ ,,,,,, 



The use of conformal symmetry is that, for a given twist, only the coefficients Cjj with 
ji + j2 < j < J — is in Eq- f!4.5p and cj in Eq. (j4.4j) for the same value of the total 
spin J can be related by EOM and/or renormalization group evolution to the leading 
logarithmic accuracy. It follows that any parametrization of DA based on a truncated 
conformal expansion is consistent with the EOM and is preserved by evolution. 

4.2 The two lowest orders 

The relevant light-cone projections of the twist-four light-cone operators corresponding 
to the three-particle pion DA are 

^11 = d+{-z)-f+igG+-{vz)u+{z) , 

V = d+{-z)^+^^gD^G,+ {vz)u+{z) . (4.10) 

In the first operator = |, ju = 1 and jg = 3/2 while in the second and the third operators 
jd = 1, ju = 1 and jg = 1. In all these cases the sum is the same, J = jd + ju + jg = 3. 
This is the minimum conformal spin corresponding to the asymptotic DA. In the fourth 
operator jd = ju = ^ but owing to the derivative jg = 2, making the minimal conformal 
spin 4. In addition, the G-parity demands that DA $|| and are antisymmetric under 
the interchange of ai ^ a2. This implies that the matrix elements of the operators with 
minimal conformal spin vanish identically in these two cases and the conformal expansion 
for <l>|| and starts one unit of spin higher at J = 4 and J = 5, respectively. 



20 




Figure 3: The two-particle twist-four DA (f)\ (left) and —(P2 (right) based on the renormalon 
model (full line) and model of Ref. (crosses) which uses the first two orders in the conformal 
expansion with the sum-rule estimate e = 0.5. The normalization is chosen to be the same 6=1. 

The model developed in jHllISI corresponds to taking into account contributions of the 
lowest two conformal spins J = 3 and J = 4. To this accuracy 

$_L(a;i, ^2, as) = 10 5^ (en - a2)al [1 + 6 e (1 - 20:3)] , 

$I|(q;i, a2, cts) = 120e(5^aia2a3(ai — ^2) , 

^±(ai, 02, as) = 105^a^(l - as) [1 + 6e(l - 2as)] , 

\l/||(ai, a2, as) = — 405^aia2as [1 + 3e(l — 3as)] , 

S^(ai,a2,as) = 0. (4.11) 

The terms ~ 5^ are J = 3 and define the asymptotic DA while contributions ~ 
correspond to J = 4. For each spin J = 3 and J = 4 there exists only one independent 
non-perturbative parameter and from conformal symmetry it follows that they have to 
have an autonomous scale dependence. This can be checked by a direct computation. 

We have defined an alternative, renormalon model of higher-twist DA in Sect. 3 by 
requiring that it reproduces the correct normalization in Eq. ()3.1|) . It is now a matter of a 
simple algebra to expand the renormalon model in Eq. ()3.4|) in contributions of increasing 
conformal spin (j4.5p and compare with Eq. (j4.1ip . We find that the structure in Eq. (j4.1ip 
is reproduced and the parameter e proves to be independent on the choice of the leading- 
twist pion DA. One obtains 

eiRen = 7/12 ^ 0.58 (4.12) 
which is comfortably close to the QCD sum-rule estimate 

e|gj^ ~ 0.5±0.1 [6]. (4.13) 

It is important that the conformal expansion of all four DA in Eq. (j4.5|) yields the same 
value of e, which illustrates that the renormalon model is consistent with the EOM. Note 
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that 02, eta) = to this accuracy; as mentioned above, its conformal expansion 

starts with J = 5. 

The two-particle twist-four DA of Eq. ()2.3jl are related to the three-particle ones 
through the EOM, Eq. (j2.12|l . Thus, when truncated to order J = 4 in the conformal 
expansion, the renormalon model essentially coincides with the model of (HI , with the sole 
difference being the replacement of the sum-rule estimate for e in Eq. ()4.13p . by that in 
Eq. ()4.12|) . However, in the renormalon model there is a priori no need to truncate the 
expansion at J = 4: Eq. ()3.8|) represents the sum of all conformal spins. Figure E] 
compares the two-particle DA between the model of jB] and that of Eq. ()3.8|1 . In general 
the two models are close, especially for 02'*\ We note, however, a qualitative difference in 
the end-point behavior. This is the subject of the next section. 



4.3 Higher orders and the end— point behavior 
4.3.1 Three— particle distributions 

Let us start with the three-particle functions. The difference between the models in 
Eq. ()3.4|) and Eq. ()4.11|) corresponds to the contribution of higher conformal spins J > 5. 
Most striking is that the end-point behavior of the renormalon model expressions for 
small gluon momenta as = ^ is $ ~ const for all five three-particle distributions 
in question whereas for each order in the conformal expansion $||, ^I/y ~ ^3, ~ af 

and ~ 

The difference indicates that the conformal expansion is not converging. Indeed, as- 
suming the asymptotic leading-twist pion DA in Eq. ()3.6|) we derive** the following formal 
expansion of the DA in Eq. ()3.7|1 in contributions of the conformal spin J: 



$iifai 



^2 ) ai + a2 ( _2 



^ (2J-1)(J+1) ^(2,2),, . 



J=3 
00 J— I 



J- 1 

i + (-i)n (2j-i)(2j-i)(j+j-i) 



x(l - as) 



J=3 j=2 



-2 p{l,2i-l)r. 



2^3] Pi]:^ 



(J-l)J(J-l) 



a2 — Oil 

1 - "3 



(4.14) 



Taken literally the expansions in Eq. ()4.14|1 are badly divergent for any fixed value of the 
gluon momentum fraction 03 and have to be understood as distributions in mathematical 
sense, i.e. they have to be convoluted with a suitable test function. Note that for the 
"transverse" distributions the double sum disappears since the only non-vanishing con- 
tributions come from j = 3/2 (n = 0) terms, c/. Eq. ()4.5|) . Owing to G-parity the DA 
\E'^,\E'll are symmetric and are antisymmetric under replacement ai ^ a2- This 

symmetry is realized differently for the transverse and the longitudinal components: for 



**To project the different DA onto the basis one uses the coefficients of Eq. (|4.9|l . The integration over ai 
and a2 separates into a product of two independent integrals upon changing variables to r = 2(q;i + (22) — 1 
and a = {0,2 ~ Oii)/{ai + 0:2) with both integrals ranging from —1 to 1. 
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^± it is associated with the exphcit overall factor ai ±^2 while the expansion itself is 
symmetric, depending only on ai + a2 = ^~<^3- For '^'||; ^||; on the other hand, the proper 
symmetry is obtained by the selection of odd/even values of j and taking into account 

thati^^i? [x] = i-iy pj]:^ [-X]. 

Absence of convergence may be an artifact of the single-dressed-gluon approximation. 
As well known |T2|, this accuracy is sufficient to identify the position of the singularity 
w = 1 but does not distinguish between singularities of different strength so that all 
renormalons appear as simple poles^^ In the full theory they will be converted to branch 
points and instead of a pole at w = 1 for a given J there will be a sum of terms with 
different singular behavior 

1 ^ rjh 



h: 



where f5o is defined in Eq. ()2.19j] and 70'^''*'' are the eigenvalues of the leading-order anoma- 
lous dimension matrix 7q"'^ for operators 0,j{^^) with conformal spin J: 

= -1^'' OJ{^?\ 7^^) = 7j'^«./vr + . . . . (4.16) 

For large spins the anomalous dimensions are dominated by soft-gluon emission and for 
quark-ant iquark-gluon operators one expects 38j 

In J < 7^^'''^ < A/-, In J , (4.17) 

where the prefactors are nothing but color charges corresponding to the possible classical 
geometries of color fiow. The logarithmic rise of anomalous dimensions translates to the 
suppression of contributions of higher conformal spin operators at large scales 



const-ln J 



^— const-m In 



in the same way as higher-spin contributions get suppressed for the leading-twist DA [HE]- 
This suppression will improve the expansion in Eq. (j4.14p and make it convergent at very 
large scales. In this sense, the renormalon model in Eq. ()3.4|) . Eq. ()3.7p can be regarded 
as representing a worst-case scenario for the convergence of the conformal expansion. 



4.3.2 Two— particle distributions 

The large higher-spin contributions to three-particle pion DA in the soft-gluon region 
03 ^ do not necessarily yield large corrections to physical observables because the gluon 
momentum fraction is always integrated over. Two-particle DA are more directly relevant. 

^^This corresponds to picking up the leading quadratic divergence of the twist-four operators and 
ignoring possible logarithmic enhancement /suppression. Also note that we do not show in Eq. H4.15|l 
the term ~ Pi/Pq which usually appears in a sum with Jq'^'^^ /Po because we tacitly imply using the 
scheme-invariant Borel transform where u is conjugate to ln(y^^/A^) and not to l/(as(/Lt)/3o/7i'). 
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Figure 4: The first few orders (partial sums) in tlie conformal expansion of the renormalon 
model for the two-particle twist-four DA 0^^^ (left), —(j)^^ (middle) and (j)^^^ +02^^ (right). The 
thick black line represents Eq. (|3.8)) . 



We find that 



and 



4^)(n) 



in Eq. ()3.8|) in the renormalon model have the asymptotic 



,2/ 



and 



u 



u 



l^lnwCl 



u 



respectively, 



behavior in the end-point regions 

which is to be compared with the ~ m^(1 — m)^ behavior of the leading conformal spin 
contributions (asymptotic DA) in both cases [H]. 

The expansions we obtained for the three-particle DA can readily be inserted into 
Eq. ()2.12|) to yield the conformal expansions of the two-particle twist-four amplitudes. In 
the case of ^g^'' the integration can easily be performed with the result 



u\l-u) 



E 

J=3,5,7,. 



j(j_l)2(j_2)^^-3 



Away from the end-points one can use the asymptotic expansion jSH] 

pK/3).„„. m _ cos{[n +{a + (3+ l)/2]e - (2a + l)7r/4} 



cos( 



+ 0{n 



-3/2N 



(4.18) 



(4.19) 



™ [sin(e/2)]"+i/2[cos(e/2)]/3+i/2 

to see that Eq. ()4.18|1 is convergent. It is not converging uniformly at the end points, how- 
ever, which explains the logarithmic enhancement compared to the asymptotic DA (and 
the model of [HI). As noted in jB], the second derivative {d'^ / du'^)(j)^\u) corresponds to 
the non-local light-cone operator d{—z)'~f^'~fr^u{z) = d_{—z)'~f^'~f^u_{z) with both quarks 
having spin projection s = — 1/2 alias jj = ju = 1/2. It follows that the conformal expan- 
sion of {d'^ /du'^)(j)^\u) goes over Legendre polynomials, or pj°!°^(2'u — 1). This result is 
consistent with Eq. (14.1^ 

^ .o„r.,u„ (J-2)(J-l)Pj°_:?(2M-l). (4.20) 
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du"^ 



u 



'l-ufPf^^{2u-l 



For 



a similar, closed-form all-order expression is not available. However, it is 



straightforward to compute the conformal expansion order by order in J, c/. Eq. ()1.2|) : 



5 

2 L 



2 -2 
U U 



j=3 ^ 24 



uu (13 uu + 2) 



ttin notations of {dP / du'^)(t>^^\u) 



{d/du)g2{u). 
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Figure 5: Subsequent approximations by the conformal expansion truncated at spin Jmax = 
3,4,... for the integrals Iq (left panel) and Ii (right panel) defined in Eq. 1)4. 22() for cj)^^^ (n) 
(diamonds), (p'^\u) (circles) and (j)^i\u) + (l)'^\u) (crosses). The results are all normalized to 
the exact value obtained using Eq. (|.S.8|) . Note the different scales on the vertical axis. 

+ 2 {6u^ + 3u + 1) \n{u) + 2 {6u'^ + 3u + 1) uHn{u) 



J=4 



+ 



1 21 

-uu (140 + 243 MM + 42) + — (6 + 3 M + 1) \n(u) 
L40 20 



+ ^(6u2 + 3u + l)u^ \n{u) 



J=5 



+ 



(4.21) 



Note that each term in the conformal expansion is of the order of u^(l — u)^ near the 



uil 



u 



behavior 



end-points but the expansion is not converging uniformly so that the 
emerges in the sum of all spins. 

The first few orders (partial sums) in the conformal expansion of the renormalon model 



4« 



and 



are compared to the full result in Fig. |3] for the three functions: 
The convergence is worst for the latter because of partial cancellation of the leading terms. 
Absence of uniform convergence at the end points means that the limit m — 0, 1 and the 
summation over conformal spins cannot be interchanged; the partial sums Y. J^^ diverge 



as 



<P^\u) ~ (5V2)«2(i _«)V1, and <p'^\u) ~ -25^u\l 



(4), 



U 



In Jmax, respectively. 



Last but not least, we show in Fig. the subsequent approximations by the conformal 
expansion truncated at spin Jmax = 3,4, . . . for the typical integrals that one encounters 
in the description of hard exclusive processes in QCD: 

du 



Jo u 



h= i'-Hl/u)<pf\u). 
u 



(4.22) 
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The difference between the renormalon model and the model of jB] (Jmax = 4) is of order 
10-30% for the first integral and somewhat larger for the second one. 

5 Renormalon model for twist— four DA of the rho 

A useful feature of the renormalon approach is its universality: with minor modifications 
it can be applied to DA of vector mesons as well. For definiteness we consider here 
p+ meson. One difference to the pion case is that because of spin the number of DA 
proliferates significantly. We will conform to the definitions and notations of Ref. ^2] 
and in particular distinguish between chiral-even and chiral-odd DA corresponding to the 
operators with odd/even number of 7-matrices, respectively, between the quark fields. 
A second difference is that because of a sizable p-meson mass the twist-four p-meson 
DA receive the Wandzura-Wilczek-type contributions of the operators with geometric 
twist -two given in terms of the leading-twist p-meson DA with the same (longitudinal 
or transverse) polarization, and the operators of geometric twist-three that are expressed 
in terms of twist-three DA with the opposite polarization ^ The Wandzura-Wilczek 
contributions to the vector-meson DA have been calculated in Refs. |1H I12j. They have 
to be added to the "genuine" twist-four contributions considered here. In this section 
we collect the necessary definitions and summarize the results; see Appendix C for more 
details. 

In what follows stands for the p-meson momentum, = m^, and e[j^^ is the 
polarization vector e ■ P = 0. We use the notation 

ml 

p^, = Pf,- z^, P = 0, 

2p ■ z 

and define the transverse polarization vector by 

We will also use the projector onto the directions orthogonal to p and z: 

9i,v = 9f^u {p^,z^ + PuZf,) . 

^ pz 

5.1 Chiral— even distribution amplitudes 
5.1.1 Definitions 

We start by quoting the necessary definitions from Ref. [T2| and in this section consider 
matrix elements involving an odd number of 7 matrices, which we refer to as chiral-even 

^The mismatch is due to different twist definition: "dimension minus spin projection" (coUinear twist) 
for DA vs. "dimension minus spin" (geometric twist) for operators, see ^JEl for details. 
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in what follows. For the vector operator the light-cone expansion to twist-four accuracy 
reads: 



(0| di-xMx) lp+(P,A)) = fpm, 

L J fij. 







2 2 
Tft X 

H — - — mu] 



4 



1 e^-^^x 



Jo 



{5.1] 



We do not consider the axial-vector operator because its light-cone expansion only includes 
contributions of twist three, five, etc., that are not relevant in the present context. For 
brevity, in this section we do not show gauge factors between the quark and the antiquark 
fields; we also use the short-hand notation 

e = 2m - 1. 

The vector decay constant fp is defined, as usual, as 

(0|J(0)7,«(0)|p+(P, A)) = /,m,e(^) . (5.2) 

The expansion in ()5.1|) involves three Lorentz invariant amplitudes which we have to 
interpret in terms of meson DA. Definitions of the latter involve non-local operators at 
strictly light-like separations and can most conveniently be written for longitudinal and 
transverse meson polarizations separately. Following fUE], we define chiral-even two- 
particle DA of the p"*" meson as 



d{-z)^l,u{z) |p+(P, A)) = fpTTlp 



'z 



PZ Jo 



+ e 



(A) 



2 ''{pz) 



due *^^''0||(m,p|) 



due "^^^""gs^u, fip) 



(5.3) 



The distribution amplitude (f)\\ is of twist two, g^^^ of twist three and g^ of twist four. All 



three functions (f) = {(f)\\, g^f^ , gs} are normalized as 



du(j){u) = 1. 



(5.4) 



This can be checked by comparing the two sides of the defining equations in the limit 
Zp, ^ and using the EOM. 

Comparing (|5.3|) with the light-cone expansion in (|5.1|) one finds |12j 



C{u) = g,iu) + <P\\iu) - 2g^;:\u) . 



(5.5) 
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The remaining invariant amplitude A(m) accounts for the transverse momentum distribu- 
tion in the valence component of the wave function. We end up with two two-particle 
twist -four DA of the longitudinally polarized p-meson, A{u) and (73 (m) which are coun- 
terparts of the pion DA 0i and 02 (a precise correspondence will be given below). 

Three-particle chiral-even distributions are rather numerous and can be defined by 
the following matrix elements: 

{0\d{-z)gG^,^{vz)-fal5u{z)\p'^{PA)) = fpmpPa[Pue''ll - p^,e''ll]A{v,pz) 

+ fp^'l [P^^9au - Puga^2'^iv,pz) 

pz 

+ fpmpj—^Pa[Pf^Zu - PuZf,]'^{v,pz), (5.6) 
{0\d{-z)gGf,^{vz)i'yau{z)\p'^{P,X)) = fpmpPa[pue^ll-Pf,e^ll]V{v,pz) 

pz 

+ fp^lj—^PcAPt^Zy - PvZp\^!{v,pz), (5.7) 

where 

^(.,p.)^/p..e-<™)^(„,), ,5.8) 

etc., and at is the set of three momentum fractions {ai,a2,0i^} = a^, a^}. The 
integration measure is defined in Eq. ()2.11|) . 

Similarly to the pion case, all higher-twist two-particle DA of the p-meson do not 
present genuine independent degrees of freedom but can be expressed in terms of three- 
particle DA. The corresponding relations ^2] are given in Eq. ()C.1|) below. 

In addition, we introduce a new twist-four DA 



(0 \d{-z)-fpgD''Ga(3ivz)u{z) \ p+{P, \)) = fpmlp^pp '^p{v,pz) , (5.9) 

pz 



which can be viewed, equivalently, either as a three-particle quark-antiquark-gluon DA, 
or as a special case of a four-quark DA with the quark- ant iquark pair in a color-octet 
state and at the same space point, c/. Eq. ()2.16|) . 

5.1.2 Renormalon model and comparison with Ref. [12j 

The renormalon model for twist-four DA of the longitudinally polarized p-meson can 
most simply be derived from the UV-renormalon ambiguity in twist-four operators. The 
calculation is very similar to that in the pion case. One obtains in the single-dressed-gluon 
approximation 

5^y\^z^d{-z)Yl^gGp^{vz)u{z)^ = 2ck^ da{l + a) d{-y)^u{z) + d{-z)^u{y) 
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^vv^z^d{—z)^'y5gG^p{vz)u{z)j = 2cA'^ZpJ daa d{—y)^u{z) + d{—z)^u{ 
Svv^z^ d{-z)YigGf,^{vz)u{z)^ = -2cA'^ [ da{l - a) d{-y)^u{z) - d{-z)^ 

5^^\^z^d{-z)^igGpp{vz)u{z)^ = 2cA^ z^ 
6yy\z'^d{-z)gD"'Go,i3{vz)^u{z)\ = -4icA^ / daa 



da a 



d 



d{—y)^u{z) — d{—z)f:u[ 
1 



daVl + V 
1 



d{-y)h{z) 



+ 



l-v 



d{—z)^u{ 



(5.10) 



where c is given in Eq. ()2.33|) and the variables y and y are defined below Eq. ()2.38|) . To 
obtain the renormalon model, we take the matrix elements of the operators in Eq. ()5.10|) 
between the vacuum and a p-meson state, project onto the relevant Lorentz structures 
and fix the normalization by the local matrix element 



(0|J(0)(7G'^,(0)7"75^(0)|p+(P,A)) = /X4-'^C4, C4 = 0.15 ±0.10 m 



making the substitution 



cA2 ^ ^Uml. 



(5.11) 
(5.12) 



The resulting twist-four DA, expressed in terms of the twist-two one, read: 
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^(a^) 



rC4 



$(a») = -U 



(l-ai)2 (1-^2)2 

a20||(ai) _ «10||(Q2) 

L(l-ai)2 {l-a2Y\ 
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1 — «! 1 — a2 

0ll(«2 



6"" 



1 — «! 1 — 02 

^2 

1 — ai 



011 



1 - 0:2 



011 («2) 



(5.13) 



where we used the symmetry 0||(u) = 0||(1 — u). Note the similarity of this result to the 
pion case, Eq. 1)3.4^ : upon replacing in the latter 5^ by (^4 one recovers the former with 
— > <^f, $1 — > m\ — > — > and — ^ Hp. 
The two-particle twist-four DA can now be restored using EOM ^21- The calculation 
(see Appendix Kill) gives: 









/ dv 0| 


{v){ 




Jo 




V 



u 



+ + iu — v) In 



u — V 



29 



+e{u <v) — 



u + + iy — u)\n. 



V — u 











Ill") 


6{u > V 
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(5.14) 
2) and g-i to 



Continuing the comparison with the pion case K{u) is similar to 16 (0i 
2(fct)2/du^. 

Eqs. fj5.13p and ()5.14p are vahd for an arbitrary leading-twist DA. Choosing the asymp 
totic expression = 6m(1 — u) yields a simple model 



1 — «i V — (X2 

= U\a2 - ai], 
Sp(ai) = 




(5.15) 



and 



A('u) = 16C4|m 'uln('u) — Li2('u) 
^73(m) = 4C4{ln(M) + ln(M)+2}. 



u \\i{u) — Li2('u) 



2uu + 



vr 



This model should be compared with that of Ball and Braun (BB) in Ref. 
the two first orders in the conformal expansion: 

63 



(5.16) 
based on 



,BB 



40 C4 o.\a-2(y.j, 



1 + -^f (3«3 - 1) 



\1/ (oj) = 630^4^014 aia2tt3(a2 — 

63 



$^^(a,: 



sf(«0 



10C4«3(l-«3)[l + Y^^(2a3-1) 



10 C4 03(02 — tti) 
0, 



63 4 , 

1 + y^f (2^3 - 1) 



(5.17) 



where the new parameter uj^ is defined by the matrix element of the J = 4 local operator 
in Eq. (4.21) in Ref. ^21- From QCD sum rules ^2] one obtains a crude estimate: 



UJ 



4 ISR 



0.8 ±0.8. 



(5.18) 



To the same J = 4 accuracy 
200 



A^^(u) = —C,av?u^ - 42C4wf |uu(2 + Viuu) + 2m^(1 + 3m + 6m^) Inw 



30 




Figure 6: The two-particle twist-four DA of the p in the chiral-even sector, K[u) (left) and 
g^{u) (right) based on the renormalon model of Eq. ()5.16|) (full line) and model of Ref. jl2j . 
given by Eq. ()5.19|) (crosses), which uses the first two orders in the conformal expansion with 
the sum-rule estimate = 0.8. The normalization is chosen to be the same ^4 = 0.15. 

+ 2#(l + 3u + 6u2)lnM|, 
2(1 

gT{n) = --Ul-Quu). (5.19) 

To avoid misunderstanding, note that for this comparison we suppressed the Wandzura- 
Wilczek contributions of twist-two and twist-three operators to the coefficients in the 
conformal expansion in ^| and only retained the genuine twist -four contributions. 

Performing the conformal expansion of the renormalon model, Eq. ()5.15p . up to J = 4 
we recover the structure of Eq. ()5.17|) predicting: 



^4lKen = (5-20) 



which can be contrasted with the sum-rule result in Eq. ()5.18|) . Similar to the pion case, 
this number is not sensitive to the shape of the leading-twist DA. 

Figure ini compares the renormalon model with that of Ref. jT^. For (73 (m) the main 
difference is in the asymptotic end-point behavior: it is logarithmic in the renormalon 
model and constant in the model of Ref. ^2|- For A(m) the difference is more pronounced 
— the asymptotic behavior is linear and quadratic in the two models, respectively — and, 
moreover, it extends to the central region because of the very different J = 4 contributions"'" 
which are determined by Eq. ()5.20|) and the central value of Eq. fj5.18|) . respectively. 
Figure [7| makes a similar comparison but this time adding the Wandzura-Wilczek terms 
in both models. 



■I-Note that there is no J = 4 contribution for gziu). 
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Figure 7: The same as Fig. ^ but including (in both models) the Wandzura-Wilczek contribu- 
tions as given by Eqs. (4.22) and (4.23) in Ref. il2^ with the parameters of Table 2 there, except 
that we consistently use the asymptotic leading-twist DA, so a| = 0. 



5.2 Chiral— odd distribution amplitudes 
5.2.1 Definitions 

For the chiral-odd operator involving the cr^i,-matrix the light-cone expansion to twist- 
four accuracy reads [12j: 

(0| d{-x)a,,u{x) y{P, A)) = 2/J (eWp. - e^P^) /rf«e-^«^^ <P^{u) + ^Kt{u) 



1 



-m. 



(Px)2 P 

2 ^1 



Px 



^ due-*«^^CT(u)|, (5.21) 



where the tensor coupling /J is given by 



\u{ 



{0)a,MO)\p-{P, A)) = ./J(eWp. - e^P,) 



(5.22) 



In comparison, the corresponding light-cone DA are defined as 



di-z)a,M^)] Jp\P,X))=^fJ\ie^ly-e^llp,) / due-'^^^M^, pi) 



(pz)2 " 

2 /"I 



eS^^)— / c?Me-^«P^/i3(M,/i|)|, (5.23) 
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where (f)±{u) is the leading-twist DA of the transversely polarized p-meson, h\^^\u) is of 
twist -three and of no interest for our present purposes, and h^{u) is of twist-four. All 
three functions (p = \ ^3} are normalized as du(f){u) = 1. Matching ()5.23|) with 

the light-cone expansion in Eq. ()5.21|) one obtains 



Ct{u) = h-i{u) - 



(5.24) 



The remaining invariant amplitude At{u) describes the transverse momentum distribution 
in the leading-twist component of the wave function. We end up with two two-particle 
twist -four DA of the transversely polarized p-meson, At{u) and h^lu) which are counter- 
parts of the distributions A(u) and gsi^u) for the longitudinally polarized p-meson and of 
the pion DA (f)^^\u) and <^r(«) defined in Eq. (EH). The precise correspondence will be 
given below. 

The three-particle DA are even more numerous than in the chiral-even case and can 
be defined as [T^ : 



{0\d{-z)a^f,gG^,{vz)u{z)\p+{P,X)) 



pT 2 



(A) ± (A) ± (A) ± I (A) ± ■ 



(A) ± (A) _L (A) ± (A) _L ■ 

Pti^±a9pu Pfi^±p9au Pv^±a9p^i + Py^^fjQa^i 
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pz 



(A) (A) (A) , (A) 

PaPi^eYi^Zy - ppp^e^Zu - PaPveYpZ^ + ppPue'^z^ 



pz 



(A) (A) (A) , (A) 

PaPi^GllZp - PpPf^e'^Za - PaPue'^Zfj + PpPyCy^Z^ 



T{v,pz) 

Ti{v,pz) 
T2{v,pz) 

Ts{v,pz) 

T,{v,pz), (5.25) 



{0\d{-z)gG,,{vz)u{z)\p+{P,X)) = tfjml 



(A) (A) 



{0\di-z)igG,,{vz)^,uiz)\p-^iP,>^)) = ^/pV 



(A) (A) 



S{v,pz), 
S{y,pz). 



(5.26) 



Of these seven amplitudes T is of twist three and the other six of twist four; higher-twist 
terms are suppressed. 

We also introduce one more twist-four DA as follows 



i<d\d{-z)a^,gD'^Go.p{vz)u{z)\p-^{P, A)) = ifpin] 



(A) (A) 



Pp'^Uv.pz). (5.27) 



As in the other cases the twist-four two-particle DA do not present independent de- 
grees of freedom and can be expressed in terms of three-particle DA using EOM ^21; see 
Eq. (insi) below. 
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5.2.2 Renormalon model and comparison with Ref. [12j 



Computing the first UV-renormalon contribution to the operator d{—z)aai3G^u{vz)u{z) 
and taking the relevant Lorentz projections we find 

6^J d{-z)e^iy a.pgG^-'ivzMz) 



-rf e'-llz, I da{l + 2a) d{-y)(yP'u{z) - d{-z)aP'u{ 



z)e^iya''^gG^,{vz)u{z) 



icA'e^'lz, 



da (1 - 2a) d{-y)a^'u{z) - d{-z)aP'u{ 



.[d{-z)e^llz^z^a''f gG''''{vz)u{z) 



-2icA^ z'^ZpC^^^l I da a 



,[d{-z)e^llz^z,a''^ gG^''{vz)u{z) 



d{~y)aP'u{z) - d{-z)aP'u{ 



d{-y)aP'u{z) - d{-z)a'''u{ 



(5.28) 



where y and y are defined below Eq. ()2.38|1 . In a similar manner, for the operators 

d{—z)Gai3{vz)u{z) and d{—z)'j5iGai3{vz)u{z) we obtain 

rl 



5vv\d{-z)e]^^z'' gGpuivz)u{z)> = cA z'^eYfi / da d{-y)(7ai3u{z) + d{-z)aapu{y) 



5v^i\d{-z)e'^lz^'^^igG^y{vz)u{z)^ = -cA^ z"eY^^ j da Y'{-y)ao,pu{z) + d{-z)ao,pu{ 



,2 ^a.W 



(5.29) 

and finally for the operator d{—z)a^ygD"Ga(3{yz)u{z) we get 
5v^^z^d{-z)gD°'Gap{vz)a^uu{z) 

= -AicK^ I daa-^ ^ ^ d{-y)cra,,u{z) + ^ d{-z)aauu{y) .(5.30) 
Jq dall + V 1 — f J 

Taking the matrix elements of the operators in Eq. ()5.28|1 through Eq. ()5.3()|1 between 
the vacuum and the p-meson state we extract the renormalon ambiguity for these twist- 
four DA in terms of 0_l(m). Going over from the renormalon ambiguity to a model, one has 
to take into account that in the present case there exist two independent local operators 
of the lowest dimension that have proper quantum numbers: 

mO)9G,MumpHP, A)) = ^f^mieM^-^P^ - eWp,), 
(0|J(0)^7G,.(0)z75n(0)|p+(P, A)) = if,m'ff{efP, - e^P,). (5.31) 
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The parameters (J ± (J renormalize multiplicatively with different anomalous dimensions 
and from the QCD sum rules one finds [101 



Q - Ci 

cr+2 



0.2 ±0.1 
0. 



(5.32) 



The vanishing (or smallness) of the second number in Eq. ()5.32|) appears as a consequence 
of vanishing of the leading contribution to the corresponding correlation function, see 
Appendix C in [T^ . 

By comparison of the expressions in Eq. ()5.29|) we observe that the leading UV- 
renormalon contribution to the operator d{—z) [gG fj^i,+ gG fj^iyi'y^Kv z)u{z) and thus to Cj+Cj 
vanishes as well. Therefore, within the renormalon model (J = —Cj and similarly to the 
pion and the chiral-even p-meson cases the model has only one parameter. 

Collecting everything and making the substitution 



2 

,4 



(5.33) 



we obtain 



Ti{ai) 
T2{ai) 
S{ai) 



a20±(ai) _ ai0±(a2) 
(l-ai)2 

0±(ai) 



S{ai 



-2cr 



2^4 



02 
1 — «! 



l-ai) 



;i-«2)2 

_ 0±(a2) 

' (l-a2) 
'±(a;2) 



1 — ai 
L(ai) - 



1-02 



1-^2 



"2 



(5.34) 



where we used the symmetry 0_L(n) = 0_l(1 — u). The correspondence with the pion DA is 
as follows: upon replacing 5^ by 3Cj, — ^ Tf , — > -T!^, m\ — > and — > Sj. 
There is no analog for '^J^. 

Finally, the two-particle twist-four DA are restored using EOM (see Appendix IC.1|1 : 



Kt{u) 



dv 4>±{v) 



[u > V) — \ u + [u 



+6iu < v) 



h{u) = 0. 



v) In ■ 



u + (v — u)ln 



u — V 

V 

V — u 



V 



(5.35) 



Note that there is no renormalon ambiguity in hs{u) at the level of a single dressed gluon, 
and therefore in our model this DA vanishes. Also note that apart from the different 
overall normalization (5^ — ^ 3^4") At{u) is the same as (m) in the pion case. 
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Choosing the asymptotic leading-twist DA 4>±{u) = 6u{l — u) we obtain simple ex- 
pressions 



and 





= -T3(ai) 


= QC,Jaia2 


i 


1 — tti 1 


m ( \ 


r-r-\ / \ 


= 3Cj[a2 - 






= -S{ai) 






) = o 






48crj 


u ln(u) — 


Li2{u) +u 


ln(u) - Li2 



(5.36) 



— uu + 



TX 



hs{u) =0. 



(5.37) 



These results can be compared with the model of [12] . The structure of the conformal 
expansion to spin J = 4 accuracy is more complicated in this CclSG clS it involves three 
parameters. Following we write 



Ti{ai) = 120tio(a2 - ai)a2aia3, 



T2{ai) = -30al{a2 - ai] 



sqo + 2 -510 (Saa - 3) + SqiOs 



Tsioi) = -120tio(a2 - ai)a2aia3, 

1 



Ti^ai) = 3003(02 — a 



soo + 2 -510 (5a3 - 3) + soioa 



(5.38) 



S{ai) = 30al 



S{ai) = 30^3 



soo (l-tts) + sio {asil-a^) - -{al + aj) )■ + Sqi (03(1-03) - 602O1} 
soo (l-as) + S10 <jo3(l-03) -^{al + al) } + Sqi {03(1-03) - 602O1} 



where, assuming that Cj = —Cly one obtains 



Soo 



"Soo — <,4 



(5.39) 



and the remaining six coefficients involve three parameters {{Q^^^)), {{Q^^^)) and {{Q^^^)) 
defined by reduced matrix elements of local operators specified in Eq. (5.20) in |[12j: 

sio = |((g(^))) + ^((g(^))) + y((g(^))), 

Sio = -|((g«))-^((Q(^))) + y((Q(^))), 

soi = +^((Q«))-^((g(^))) + ^((Q(^))), 
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ho= + (5.40) 

Here we only show genuine twist-four contributions and suppress the Wandzura-Wilczek 
terms. To the same accuracy 

A^^iu) = UOClu'u' - ((Q(i))) + ^ ((Q(3)))) [uu{2 + 13uu) 

+ 2m^(10 - 15u + Qu^) \nu + 2u^(10 - 15m + 6u^) \nu , 

/i^^(m) = 0. (5.41) 

The particular model suggested in [T^ makes use of the QCD sum-rule estimate 

{{Q^'^))\sK = -0-15 ± 0.15 [40J , ((g(^) ))\,^ = ((g(^) ))\,^ = O . (5.42) 

On the other hand, starting with the renormalon model in Eq. ()5.36|) and Eq. ()5.37|) 
and isolating the J = 4 contribution^ we obtain, using Eq. ()5.32|) . 



{{Q^'^))\K.n = y ((^?^'^ ))lRe„> HQ'^'^ ))lRen = 'd = -CIO ± 0.05, {{Q^'^ = . 

(5.43) 

Note that the renormalon model prediction for ((Q*-^-*)) is consistent with the sum-rule 
estimate within errors and the main difference is that {{Q^^^)) is non- vanishing. With the 
numbers from Eq. (j5.43p the J = 4 contribution to At{u) is enhanced by roughly factor 
four compared to the sum-rule estimate, but is still smaller than the leading J = 3 term. 

We conclude by a numerical comparison of the two-particle DA At{u) in the renor- 
malon model, Eq. ()5.37p . and the model of |^ in Eq. 1)5.411) with QCD sum-rule estimates 
of the parameters. As in previous cases, the difference is most pronounced in the end-point 
regions where the two expressions have different asymptotic behavior. 



6 Conclusions 

In this paper we have presented for the first time a systematic analysis of twist-four meson 
distribution amplitudes that goes beyond the first few orders in the conformal expansion. 
Our analysis is based on the study of the high-order behavior of perturbation theory in the 
single-dressed-gluon approximation which is equivalent to the study of one-loop power 
divergences of the contributing twist-four operators. In general, this calculation supports 
the conjecture that the shape of higher-twist DA is not far from the asymptotic form 
which is associated with operators with the lowest conformal spin. However, we find that 



§Note that T2 = T4, implies that {{Q^^^}) = while Ti = -T3 leads to {{Q^^^}) = 4(((9^^^)) 
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Figure 8: The two-particle twist-four DA At{u) without (left) and with (right) the 
Wandzura-Wilczek terms. Each figure compares between the renormalon model (full line) 
and model of Ref. ^2] (crosses) which uses the first two orders in the conformal expansion 
with the sum-rule estimate of Eq. ()5.42|1 . Note the different scales on the vertical axis. 

the conformal expansion of the renormalon model does not converge uniformly at the end 
points. Consequently, the end-point behavior corresponding to the situation where one 
valence quark is soft, is qualitatively different between the sum over all spins and any 
fixed order in the conformal expansion. As the principal result, we obtain that the two- 
particle twist-four DA that describes the k'j_ distribution of valence quarks in the meson 
— (j)^i\u) for pion, A{u) and A-^(m) for p-meson — has the same linear falloff ~ u{l — u) 
at M — i> 0, 1 as the leading-twist DA, compared to the quadratic behavior ~ m^(1 — m)^ 
of the asymptotic DA (lowest conformal spin). Taking into account the limitations of our 
analysis this has to be understood as an upper bound. The existence of such a bound is 
important for proofs of factorization theorems. 

An attractive feature of the renormalon approach is that it allows one to construct 
simple models of higher-twist DA with minimum number of non-perturbative parameters. 
In this paper we constructed such models for the pion and for the p-meson with both the 
longitudinal and the transverse polarizations. The corresponding expressions are given in 
Sect. 15.1.21 and I5.2.2t respectively. In each of these cases the entire set of twist-four 
DA is determined in terms of the leading-twist DA with just one free parameter, the 
overall normalization, corresponding to the matrix element of a certain local operator. 
This approach presents a viable alternative to the models of Refs. jUJ [T2j based on the two 
first orders in the conformal expansion. The spread between the predictions of these two 
approaches is a fair measure of uncertainty in our present understanding of higher-twist 
effects. 

In addition to giving, for the first time, an upper bound for the possible contribution 
of the operators with large conformal spins J, the renormalon model can also be used 
to estimate the next-to-lowest spin J = 4 contributions. The corresponding estimates 
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are given in Eq. ()4.12|) . Eq. ()5.20j) and Eq. ()5.43|) for the pion, longitudinal p-meson and 
transverse p-meson, respectively. These estimates are probably more reliable than the 
corresponding QCD sum-rule results, as it is known that the sum-rule approach does not 
work well for operators containing derivatives. In particular for the longitudinal p-meson 
there is a significant difference, compare Eq. (|5.2Up and Eq. (|5.18|) . 

The renormalon approach can be applied in a straightforward manner to study yet 
higher power corrections, of twist six and above, which are otherwise inaccessible. Higher 
infrared-renormalon ambiguities in the coefficient functions of Eq. 1)2.211) at Wq = 2, 3, etc. 
corresponding to twist six, eight, etc. can be extracted from Eq. ()2.22|1 and Eq. ()2.23|) 
in full analogy with the calculation of the leading renormalon ambiguity in Sect. 12.2.21 
Assuming for simplicity the asymptotic leading-twist DA, the ambiguities in the leading- 
twist part of Gi read: 

6^=' {Ci ® 04 = (-A^A^)' 2 [u' \n{u) + ln{u) + uu] , (6.1) 

2po 4^ " ' 

15 

ST' {Ci ® 0.} = "9^^ (-A^A^)' ^uu [31 + 199UU - 2u'u'] . 

In general, we find in the renormalon approach that the asymptotic behavior of the form 
~ uu persists in the case of Gi to all twists^. For G2, on the other hand, we find at any 
twist higher than four an asymptotic behavior of the form ~ u'^u'^. This result suggests that 
the DA of all twists probably have the same universal power behavior in the end-point 
regions. This strong conjecture implies, in particular, that the twist expansion breaks 
down owing to the increasing singularity of higher-twist coefficient functions. 

The present study can be extended in several respects. From the theoretical point of 
view the renormalon calculation that uses the modified gluon propagator in Eq. ()2.20|) 
can be understood as changing the scaling dimension of the gluon field. It is, therefore, 
tempting to try to reconcile the renormalon model with conformal symmetry of QCD with 
modified conformal spin assignment for the fields. Alternatively, since the renormalon 
model predictions effectively reduce to the analysis of quadratic divergences of twist-four 
operators, they may have the symmetry of QCD in two dimensions. Finally, there is 
a general problem of going beyond the large-/3o approximation and taking anomalous 
dimensions into account. 

To summarize, we believe that the renormalon approach is useful for understanding the 
structure of higher-twist contributions in hard exclusive processes and allows one to obtain 
quantitative estimates. Phenomenological applications are numerous but go beyond the 
tasks of this work. 

^An exception is twist six (the first line in Eq. 1)6. 1(1 ') where the behavior is ~ u^u^ owing to a complete 
cancellation between diagrams 1 and 3 in figure ^ 
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Appendices 

A Single— dressed— gluon calculation of the twist— two 
coefficient function 

In this Appendix we present the detailed calculation of the diagrams in Fig. Q to arrive at 
the Borel-regularized expression for the leading-twist coefficient function of the non-local 
operator in Eq. ()2.H) . The calculation is done in the Feynman gauge and for brevity we 
do not write explicitly the Wilson line connecting operators at different points. 

We begin with Diagram 1 in Fig. ^ where the gluon is exchanged between the quark 
and the Wilson line. The latter is defined by Eq. ()2.2|1 . The calculation of this diagram 
plus its symmetric counterpart, where the gluon line is attached to the other quark, yields: 

(p2|T{J(x2)7.75m(xi)}|j9i>^ = — — ^ / dwei^{-ky / da e'^^^^^-^^^^'^ 

Po Jo Jo 

X [dp2lc7ltilulbUpA''^ii.P2, aA) + (ip27;.757/.7<T?^pi A'^I^(j9i, aA)] , (A.l) 

where dp^ and Up^ are quark spinors, /S. = X2 — Xi and the momentum integral Ii (p. A) is 
given by 

riv z) = /^'^c-^A iP + kr 

' J (27r)4 + A;)2 

1 T{-w) 



167r2r(l + w) 



(^^^ jr'rf55-e*(i-fe)pA ^^^p _ 2w^y . (A.2) 



Inserting the result of the momentum integration into Eq. ()A.1|) . contracting the indices 
and using the Dirac equation for massless quarks gives 

X "(-z6p2 ■ A + ^) e^'^^^-^)^-^ + ■ A + ^) e^'^^^-^^^'-^j (A.3) 

xe*(P^"^-^i"^)47.75Mp,. 

The terms proportional to pi(2) ■ A can be removed using —ip ■ Ae*"'-^"''-'^'^ = "^e*"^^"*''^'^ 
and then integrating by parts. Now the dependence on the external momenta is only in 
the spinors and the phase and can be absorbed in external quark states, so that the result 
takes the form: 

{P^\^^{di^^)^^lM^^)}M^ = -^J^ d j^daj^ dp 
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X J daa^'"-^^~25{a)6{/3) + J dbb"" [5{P)5 {a - a{l - b)) + 6{a)6 {(3 - a{l - b))] 

X {p2 \d{x2 + aA)-f^-f5u{xi - pA)\pi) . (A.4) 

Here we introduced for later convenience two new variables a and jS. The integration over 
a can be taken and after rearranging the terms the result can be represented as an OPE 



T{d{x2)^ul5u{xi)}-i 



2A 



' 5,„ V{-w) 



Jo 



2 a2\ /•! 



da 



l-a 



dp 



X d{x2 + aA)7''75M(xi - /?A) , 



(A.5) 



where it is understood that only leading-twist operators are retained on the right-hand 
side, 



1-1-/3 

f^'"\p) = p^""-^ / db{l-b)-^'"b 

Jo 



l + w 



2F^{[2w,w + l],[2 + w];l-(3), 



(A.6) 



and the "+" prescription is defined, as usual, by Eq. ()2.26|) . 

Next we consider Diagram 2 where the gluon is exchanged between the quarks. This 
contribution reads: 



{p2 \T{d{x2)'yul5uixi)}\pi)^^ 



Po 



where the integral l2"{pi,P2, A) is given by the expression: 



(A.7) 



ir(pi>P2, A) 



d^k 



-ikA 



(2vr) 



(pi + kY{p2 + ky 
{k'^y+'"{pi + kf{P2 + kf 



Ti-w] 



327r2r 1 + w V 4 



A^ 



2 \ /•! 



da 



d(3{l-a-(3y 



(A. 



X e 



i(op2 A+/3pi A) 



A^A'"' 



Here the dots represent terms proportional to at least one external momentum p^^}, 
which do not contribute to Eq. ()A.7|) by virtue of the Dirac equation. Inserting Eq. ()A.8|) 
in Eq. ()A.7|) . contracting the indices and absorbing the dependence on the momenta pi 
and p2 in external quark states, one immediately obtains 



T{(i(x2)7,.75M(a;i)}2 



2A 



Jo 



' I. r(-^) /-A^A- 
rfwes — 

r l + w 



2 a2\ n 



da 



l-a 



d(3 
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X d{x2 + aA)7''75u(xi - pA) . 



(A.9) 



For the last step it was important to consider a non-forward matrix element (pi 7^ P2), 
because otherwise one could not identify how the quark-field operators get shifted. 

The last contribution comes from Diagram 3 since Diagram 4, describing the self energy 
of the incoming quark, has no scale and thus vanishes. We obtain: 



T{d{x2hul5u{Xl)}3 

where 

l3(A2) = I da I db 



A' dwe-^^ i-ATUA)dix2h,iMxi) 



(A.IO) 



d^k 



-i{a-b)kA 



-W] 



327^2(l-2^^;)^(l + w) V 4 



2\ "lu-l 



(A.ll) 

Collecting the contributions in Eq. ()A.5|) . Eq. ()A.9|) and Eq. ()A.10|) we get a gauge- 
invariant result for the OPE of the T-product of quark fields to all orders in the strong 
coupling in in the large-/?o limit: 



T{d{x2)-ful5[x2, Xi]u{Xi)} 



2A 



(iwea' 



-w] 



Jo 



Til + w) 



2 \ 2\w 



-A'A 



da 



dp 



{1-a-py 



(A.12) 



l-2w 



5{a)5{p)\ X c/(x2 + aA)7''75[x2 + aA,xi -/3A]m(xi -/3A), 



where the terms in the curly brackets are grouped as they appear from individual diagrams 
in Feynman gauge: the first line corresponds to Diagram 1 (vertex correction) in Fig. 
and its symmetric counterpart, the second line to Diagram 2 (box diagram), and the third 
line to Diagram 3 (self-energy like correction to the Wilson line). Disentangling the two 
Lorentz structures we end up with the answers for the unrenormalized coefficient functions 
given in Eq. (IT^ and Eq. (IT^ . 



B The operator dj^^gDaG^^u: UV divergence and 
EOM relations 

The calculation of the UV-renormalon ambiguity of the operator in Eq. ()2.16p goes along 
the same lines as in Sect. 2.3 and is fully analogous to the similar calculation (with no 75) 
in [T71 131j. For the matrix element between off-shell quark states (Fig. 12)) we get 

{q2 \d{-z)gD-Go,p{vz)ru{z)\q,) = ^^!_Ze-(^i+''^)- {g'^Pg^^ - g^^gl) (B.l) 
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dwel^i-A'^r 47 Vr^,, hapiq2, (1 + v)z) + 417^7%, h^^iq,, (1 - v)z) 



where the momentum integral is 

hap{.q,z) -- 



:b.2) 



Computing the integral, extracting the w = 1 residue, specifying V = 7^75 and projecting 
with x^x^ we obtain" the result in Eq. p.39|) . Taking the matrix element of Eq. p.39|l 
between the vacuum and a pion state we end up with the renormalon ambiguity of the 
DA H^(aj) given in the last line of Eq. ()2.4()j) . 

Going from an ambiguity to a model involves a replacement of the large-/3o renormalon 
residue by a physical non-perturbative parameter, a certain local matrix element. From 
[6j it is known that the normalization of all four DA $_|_ y and is controlled by a single 
non-perturbative parameter 6^ for J = 3. In the renormalon model, contributions of 
higher spins J = 4, 5, ... to $_|_ y and \&_|_ y are fixed uniquely in terms of 6"^ so that no new 
parameters appear. We are going to argue that within this construction the normalization 
of is also fixed uniquely by EOM that relate it to the J = 5 contributions to the other 
DA. Thus, all five DA are given in terms of 5^ and the leading-twist pion DA. 

The necessary constraint can be derived from the operator identity for the second 
derivative 



dxadx" 



d{—x)V u{x) = —d'^d{~x)T u{x) + gd{—x) aG{—x) V + VaG{x) 



u[x] 



2igx^ 



d 



dx 



dvv d{—x)TG u^{vx)u{x) + 2igx''d^ / dv d{—x)TG,y^{vx)u{x) 



-1 



+ 2g dv dt {1 + vt)d{-x)Tx>'x''Gf,p{-vx)G''^{-tx)u{x) 



-1 



-igx" I dv{l + v')d{-x)T[D^,G^'^]{vx)u{x) 



(B.3) 



given in Eq. (A. 9) in jT2] (see also [2^1 )• Here crG = aapG"^ and stands for the 
derivative with respect to the total translation. 

Taking the matrix element of (jB.3|l between the vacuum and a pion state and using 
the definitions of the DA we obtain 
»i 



2i 



{4:-2ipz{2u- 1)) + (-3 + ipz {2u- 1)) 



Dai e 



-ipz (ai-Q2) 



pz"^ a| 



1 ) ?,m{pz 0:3) 



cos{pz 0:3) 
pz ^3 



+ 



"3 



I ai 



pz a\ pz a1 



siYi.{pz 03) H (tti — ^2) cos{pz 03) $11 (ctj 

as 



(B.4) 



"We performed integration by parts over a to eliminate q ■ z factors in order to convert to operator 
notation. 
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^ VV"3 



pz'^ oil 



sm{pz 0:3) 



cos(]9z as) 
pz al 



where we omitted the contributions of the two-gluon operator from the right-hand side. 
These can be systematically put to zero to our accuracy as they start contributing at 
higher order in the flavor expansion. 

Expansion of Eq. ()B.4|) in powers of pz yields simple relations between integrals that 
involve all five three-particle DA. The odd powers are trivial: both the right-hand side 
and the left-hand side vanish by symmetry. The first two non-trivial relations are: 



du 



-8 0i(m) + 6 02(m) 



Dai 



(B.5) 



du (u 



u? 



-16 0i(m) + 1O02(m) 



Dai 



+ - (ai - 012) (7 (ai + a2) - 10) - - (ai - a2)(ai + ^2 - 1) $±("0 



+ (1 - 2 02 - 2 «! + 2 «^ + 2 (^11 



2^±(«.)) 



The first relation does not involve S^(Q;j) and is satisfied identically both by the model of 
[H] and the renormalon model. The second relation gives the required constraint for the 
normalization integral j Vai {ai — 0:2) ^^^(ai) in terms of the other four DA. It is easy to 
verify that in order to satisfy this constraint one must assume in the last line of Eq. ()2.40|) 
the same replacement cA^ — > 5^/6 as in the other DA. 



C Cancellation of renormalons for the p— meson am- 
plitudes 

Here we want to demonstrate cancellation of IR renormalon ambiguities in the leading- 
twist coefficient functions with the UV-renormalon ambiguities in the matrix elements of 
twist -four DA for the case of exclusive amplitudes involving a vector p-meson. Similarly 
to the pion case, we consider the simplest example: a gauge-invariant T-product of quark 
fields sandwiched between the vacuum and the meson state. 



C.l UV renormalons in two— particle DA 

To begin with, we calculate the UV-renormalon ambiguities in the two-particle DA of 
twist four. These can be obtained from the three-particle ones using EOM. 
The specific EOM relations we need in the chiral-even sector are |l2'j : 

d 1 
g^{u) = - 2— / dai / da2— [2<^{ai) + ^(a^)] 
du Jo Jq as 

pu pv nu pv p1—v 2 

A(n) = 32 dv dw {gl{w) - (p^iw)) + 32 dv dai / rfa2— [2$(ai) + ^(a^)] , 
Jq Jo Jo Jo Jo 
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+8 / dai 



l-u 



da- 



ua2 — (1 — 



[2$(a,) + ^{ai)] . 



(C.l) 



This implies that the ambiguities of (73 (m) and A{u) due to UV renormalons at w = 1 are 
given by 

Suv {gsiu)} = -2— dai daa— 5uv{2<l>(ai) + 
du Jo Jq as 



u pv 



JO 



l-v 



(5uv{A(m)} = 32 / dv I dai I rfa2— 5uv {2$(ai) + ^(oi)} 



«3 



+8 / dai 



l-u 



da- 



ua2 — (1 — u)ai 



a% 



5uv{2$(a.) + ^(«0}, (C.2) 



where the ambiguities of the three-particle DA can be read from Eq. ()5.13p using the 
inverse substitution in Eq. ()5.12|1 . For gs{u) we get 



4c 



A2 


d 


9 

Trip 


du 




d 




du 



dai 



l-u 



da-> 



4>\\{ai) 0||(a2) 



(C.3) 



{l-u) / da — — - u 

[1 — ay 







{l-af 



where we evaluated one of the a integrals. Taking the derivative in respect to u and using 
the symmetry = 0||(1 — u) yields the desired result: 



A^ r /-i r 1 1 ■ 

{g^iu)} = -Ac— { / dv(t)\\{v) e{u > v) — + e{u < v) — 
In turn, for the DA K{u) we get the following expression: 



011 (w) 



uu 



(C.4) 



A2 



5uv{A(m)} = —16c — ^<4 I dv I da 



U PV 



^0 



l-v 



dan 



011 ("1) </'||("2) 



+ / dai 



l-u 



da' 



ua2 — (1 — u)ai 
\ — ai — a2 

Taking one of the a integrals in each line yields: 



[{l-aif {l-a2f\ 

0||(«l) _ 0||(«2) 

{l-aif {l-a2f 



(C.5) 



5uv{AH} = -16c^<'4 / dv 

mz 



(1 — f ) / da 



a 



(I -a) 







da \ iu — a) In 



u — a 
1 — a 



+ u{l -u] 



hi 



a] 



{l-af 



l-u 



da I (1 — M — a) In 



1 — u — a 
1 — a 



+ u{l - u) 



{I -a) 



011 («) 
{I -a) 



(C.6) 
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Finally, integrating by parts over v in the first line and rearranging the terms we obtain: 



(5uv{A(m)} = 16c 



A2 



m 



p JO 



dv(j)\i{v)< 6{u > v) — 



u + u + (u — v)ln 



u — V 



+e{u <v) — 



u + u + [v — u) in 



(C.7) 



In the chiral-odd sector the EOM relations between the two-particle and the three- 
particle DA of twist four are given by |,12j : 



d 
I 

du 



dai 



l-u 



da-} 



S{ai) (T2(ai) - T^{ai)) 

"3 



an 



dv{2v-l)[(l>^{v) + h3{v)] + 8 dv dw [hiw) - <i>±{w)] (C.8) 



+4 / dai 



l-u 



dan 



1 ai-a2-{2u-l) 
— S{ai) 2 (-?2(«i) - Jsltti)) 



as a^ 

As a consequence, the ambiguity of h^lu) and At{u) due to UV renormalons is related to 
that of the three-particle DA S{ai), T2{ai) and Ts{ai). Using the results for the ambigu- 
ities of the three-particle DA in Eq. ()5.34j) with the inverse substitution in Eq. ()5.33|) we 
obtain 



A^ d 

Svvihsiu)} = 4c^^ / dai 



1 — u — a2 



da2 
1 



u — ai 



(1 — «! — ^2)^ 1 — «! 

0±(«2) 



0±(«l) 
(l-«l) 



_(l-ai-a;2)^ 1 - 0:2 J (1 - 02) J ' ^^'^^ 

The square bracket in the first line vanishes upon taking the 02 integral and that in the 
second line upon performing the ai integral. Thus there is no UV renormalon ambiguity 
(at w = 1): 

S^y{h,{u)} = 0. 

For At{u) on the other hand 



2 /•« 



6 



UV l^T 



u)} = 8c — ^ / da 



m 



p JO 



l-u 



dai 



^ _ 2{u - ai) _^ (w - "i)(l - «iy 



(1 - «! - ^2) (l-ai-a2)^ 



+ 



2(1 - u - 0:2) (1 - n - a;2)(l - "2) 



(1-0:1-02) (1-01-02)^ J (1 - "2)^ 

Taking one a integral we obtain 



5uv {At(m)} = 16c 



A' 



2 n 



m 



dv (j)\\{v) 



p JO 



{u > v) 



{1 — u) + {u — v) In 



(C.IO) 

0±(ai) 
(l-«i)2 

(C.ll) 

u — V 



l-v 



+6{u < v) — [ u + {v — u)\n 



V — u 



(C.12) 



46 



C.2 IR renormalons in coefficient functions 



The IR w = 1 renormalon ambiguity in the all-order perturbative calculation of the 
leading-twist coefficient functions can be obtained from a generalization of the OPE rela- 
tion in Eq. |02|l for the case of an arbitrary Dirac matrix T between the quark fields: 

nl nl—a. 

S,^,{T{d{x2)Tu{xi)}} = -cA^A^ da d(3 

Jo Jo 

X I {f^\(3)S{a) + [d{x2 + aA)ru{xi - /3A)] (C.13) 

-«-/?) (^g,, + 2^^) [d{x2 + aAh^r^Y^M^i - /9A)] }, 

where we replaced the w-integral by vr times the w = 1 residue and suppressed the gauge- 
link; c was defined in Eq. ()2.33p . 

For vector mesons there are two relevant structures, chiral-even F = 7^, and chiral-odd 
r = <Jfj,u- We will discuss these two sectors separately. 

Chiral— even amplitudes 

The matrix element of the T-product of quark operators between the and the vac- 
uum state contains three Lorentz structures which we parametrize using the structure 
function Ff. 

{0\T{d{-x)jM^)}\p+{P, X)),2 = f^m, fdue-^^^^-^-^^ (C.14) 

Jo 

The OPE to twist-four accuracy reads: 

2 2 

Fi {u, x^- /i^) = Cf^ ® (n) + A(u) , 
F2{u,x^-p^)=&^^®gt\u) + ..., 

F3(u,xV') = ®0||(u) + C(u). (C.15) 

Here is the leading-twist DA of the longitudinally polarized p-meson, g^l\u) is the 

DA of twist-three corresponding to the contribution of the transversely polarized p-meson 

and the functions A{u), C{u) represent higher-twist contributions. At leading order in as 

C^\u,v) = C2'^\u,v) = 6{u — v), alias Fi{u,x'^; p?) = and F2(m, = g^^\u), 

(2) 

while the remaining twist-two coefficient function C3 vanishes. 

In order to calculate the IR-renormalon ambiguity in the leading-twist part of the 
amplitudes in Eq. ()C.14|) we take the appropriate matrix element of Eq. ()C.13|) retaining 
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the leading terms in Eq. ()C.15|) . The result reads **: 

-16cA2 



5iR,<!cf)®0||( 



u 



mi 



•1 rl-a 

dv / da 



dp 6{v{l — a — P) + a — u) (p\\{v) 



X + + 2(1 - a - /?)] 



-16cA 

mi 



2 /•! 



(if 0|| (f ) 6'(u < f ) — I M + + (f — u) In 



6'(m > v) — i u + + (u — v)\n 



V 

u — V 

V 



2 ^2 /•! /•! 



4cA^ ^2 

— ^"T^ f dv j da I dP{l — a — P)5{v{l — a—P) + a — u)4)\\{v) 
mp du Jq Jq 

4cA2 



(if 



1 1 ■ 

6'(m > f) — + 6'(m < f) — 

V V 



uu 



. (C.16) 



Comparing these expressions with the UV-renormalon ambiguities in the twist-four con- 
tributions for K{u) and C(m) in Eq. ()C.7|) and"!"^ Eq. ()C.4|) . respectively, we observe that 
the ambiguities in the structure functions Fi cancel out, as expected. 

Chiral— odd amplitudes 

The calculation for chiral-odd amplitudes goes along similar lines. We parametrize the 
matrix element in terms of three more structure functions 



{Q\d{-x)op^u{x)\ 



To twist-four accuracy 

?T/„. Jl. ,,2^ _ r^T{2) 



Ft{u,x';fi')=Ct 
F^{u^x',^?)=Cl^^ 
Fiiu,x';f^')=C^^'^ 



2 2 
TTl X 



h[!\u) + ... 



.3 ^Mu) + CTiu), (C.18) 

where 4>_l{u) parametrize the leading-twist part. Repeating the procedure of the previous 
section we obtain the following ambiguity: 

y^2 /•! /•! /"l-a 

ml JO 



. a2 rl rl i-l—a 

5iR \ ® (t)±{u) \ = -16c— / dv da dp 5 {v{l - a - p) + a - u) M'>^) 



**The IR ambiguity of €2^^ (8) g^'' is more difficult to obtain because the twist-three matrix element 
vanishes for on-shell massless quark states. This ambiguity must be compensated by contributions of 
twist-five operators and is of no interest for our purposes. 

Taking into account the relation of Eq. (|5.5|l between C(u) and 53 (u). 
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-16c- „ 

mi 



x(/«(/?)<5(a) + /«(«)5(/5)) 

A2 ^ 







/ dv(j)T{v) 


^(n<.)l( 


Jo 





V — u 

V 



\d\U > f ) — M + (U — f ) 111 



^.^{Cl'-'-^ ®4>i-{u)\ =0, (C.19) 

The absence of an ambiguity in which is to be associated with Cr(M), resuhs 

from the fact that Diagram 2 in Fig. Q vanishes owing to the identity 7aO"^!/7° = 0. The 
absence of this ambiguity is expected based on the fact that h-^iu) has no corresponding 
UV-renormalon ambiguity — see Eq. ()C.10|) — and the relation between Ct(m) and h^iiu) 
in Eq. ()5.24|) . Comparing the result for 5ir{C^'-^'' ® </>±(ti)} with the UV-renormalon 
ambiguity in A^('u) in Eq. ()C.12|) we observe that the ambiguities in the structure function 
F'^ cancel out. This completes the calculation. 
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